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ABSTRACT 


Over  the  years,  there  has  been  significant  variation  in  the  fihnwise  steam 
condensation  data  at  NFS  on  horizontal  low-integral  finned  tubes.  With  a  view  to 
increasing  the  accuracy  of  the  data,  inserts  were  used  inside  the  tubes  to  reduce 
inside  thermal  resistance;  however,  significant  discrepancies  then  occurred  in  the 
calculated  outside  coefficient  when  compared  to  data  taken  without  an  insert.  These 
discrepancies  arose  due  to  the  data  reduction  technique  which  assumes  a  known 
inside  heat-transfer  resistance  and  subtracts  this  from  a  measured  overall  resistance. 
If  the  assumed  value  on  the  inside  is  inaccurate,  then  the  outside  value  is  equally 
inaccurate. 

The  present  work  uses  an  instrumented  smooth  tube  to  obtain  accurate  inside 
heat-transfer  correlations  both  with  and  without  inserts  and  uses  these  to  obtain 
accurate  outside  coefficients  for  a  family  of  uninstrumented  finned  tubes  with  a  view 
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I.  INTRODUCTION 


A.  BACKGROUND 

A  reduction  in  size  and  weight  of  aU  types  of  heat  exchangers  aboard  Naval 
vessels  will  allow  more  efficient  use  of  space.  The  benefits  might  include  greater 
equipment  accessibility  for  maintenance  or  greater  heat  exchanger  capacity  (without 
a  corresponding  increase  in  size  and  weight)  with  a  corresponding  increase  in  fuel 
efficiency. 

For  the  past  ten  years,  the  Naval  Postgraduate  School  in  collaboration  with  the 
David  Taylor  Research  Center  and  the  National  Science  Foundation  has  conducted 
research  that  is  directed  at  the  development  of  smaller,  more  efficient  steam 
condensers.  Improved  designs  can  result  in  significant  space  savings,  always  a 
primary  concern  on  Naval  vessels,  especially  submarines. 

Uncertainties  in  past  data  using  steam  were  apparently  due  to  the  lack  of 
detailed  information  about  the  inside  heat-transfer  correlations  used  during  the  data 
reduction  process.  Previously,  the  standard  Sieder-Tate  correlation  was  assumed  to 
be  valid  for  the  inside  heat-transfer  coefficient,  but  it  may  not  be  the  best  correlation 
to  use  with  the  particular  test  arrangement  used  in  this  research  program. 

A  large  amount  of  enhanced  condensation  data  has  been  collected  in  previous 
studies  at  the  Naval  Postgraduate  School  on  more  than  90  different  condenser  tubes 
of  varying  fin  height,  fin  spacing,  and  tube  material,  with  a  view  to  finding  an 


1 


optimum  fin  geometry  for  both  steam  and  refrigerant  condensation.  Satisfactory 
results  have  been  obtained  with  refrigerant  data  using  R-113.  However,  some 
troublesome  questions  of  possible  contamination  and  instrument  inaccuracy  still 
remain  with  the  steam  data.  It  is  felt  that  in  order  to  address  these  questions,  a 
fundamental  evaluation  of  the  heat  transfer  apparatus  on  which  this  data  was 
collected  and  of  the  data  reduction  process  was  considered  appropriate.  Current 
data,  taken  on  a  carefully  cleaned  and  calibrated  apparatus,  could  be  compared  to 
previously  recorded  data  and  a  determination  as  to  its  validity  and  reproducibility 
could  be  made.  In  addition,  a  thorough  evaluation  of  the  best  inside  heat-transfer 
correlation  would  lead  to  more  reliable  steam  condensation  results. 

B.  CONDENSATION 

Condensation  occurs  when  a  vapor  is  cooled  below  its  saturation  temperature, 
or  when  a  vapor/gas  mixture  is  cooled  below  its  dew  point.  Surface  condensation 
occurs  in  condensers  when  a  cooled  surface  (kept  at  a  temperature  below  the 
saturation  temperature  of  the  vapor)  contacts  the  vapor.  The  vapor  molecules  that 
contact  such  a  surface  stick  to  that  surface  and  condense  into  liquid  molecules. 
Condensation  may  occur  in  one  of  the  following  modes:  filmwise,  dropwise,  or  mixed 
mode  (a  combination  of  filmwise  and  dropwise)  condensation.  In  the  filmwise  mode, 
the  liquid  wets  the  cold  surface  to  form  a  continuous  film.  If  the  liquid  does  not  wet 
the  surface  but  instead  forms  discrete  drops  on  the  cold  surface,  dropwise  or  mixed 
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mode  condensation  will  occur  and  is  often  caused  by  some  form  of  contamination. 
[Ref.  1] 

The  condensate  forming  on  the  tube  surface  offers  a  resistance  to  heat  transfer 
between  the  vapor  and  the  surface,  which  increases  with  the  thickness  of  the  liquid 
layer.  Even  though  dropwise  condensation  results  in  much  larger  heat-transfer 
coefficients  than  filmwise  condensation,  it  is  difficult  to  maintain  a  stable  dropwise 
condition  over  prolonged  periods.  Therefore,  in  most  cases  condenser  design 
calculations  are  based  on  the  assumption  of  filmwise  condensation,  resulting  in  lower 
heat- transfer  coefficients  and  more  conservative  designs.  [Ref.  1] 

In  a  condenser,  the  coolant  side,  tube  waU,  and  vapor  side  thermal  resistances 
control  the  heat  transfer  rate  from  vapor  to  coolant.  Also,  for  experimental  work 
we  always  use  clean  tubes,  but  in  real  condensers  tubeside  fouling  can  play  an 
important  role  an  increasing  the  coolant  side  resistance.  The  magnitudes  of  these 
resistances  depend  on  the  fluid,  tube  geometiy,  and  flow  conditions  on  the  vapor  and 
coolant  side.  For  steam  condensation,  it  is  the  coolant  side  thermal  resistance  which 
tends  to  dominate.  Methods  to  lower  this  inside  resistance  include  the  use  of  inserts 
or  roped  tubing  to  promote  turbulence,  thereby  raising  the  convection  heat-transfer 
coefficient.  However,  such  modifications  lead  to  increased  pressure  drop  through 
the  tubes,  which  must  be  compensated  for  by  providing  extra  pumping  capacity. 
Heat  transfer  through  the  tube  wall  is  conductive  and  is  fixed  once  tube  thickness 
and  material  are  selected.  The  vapor  side  resistance  is  due  to  the  condensate  film 
which  forms  on  the  outside  of  the  tube.  For  filmwise  condensation,  the  outside 
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resistance  can  be  reduced  by  the  addition  of  low  integral  fins.  These  have  the  effect 
of  not  only  increasing  the  outside  surface  area  of  the  tube,  but  also  of  thinning  the 
condensate  film  around  the  fins  due  to  surface  tension  forces.  Too  small  a  fin 
spacing  may  result  in  condensate  flooding,  whereas  too  large  a  fin  spacing 
approaches  the  smooth  tube  case;  there  should  be  an  optimum  fin  spacing 
somewhere  in  between  these  two  extremes.  Horizontal  fin  spacing  is  therefore  of 
prime  importance,  and  finding  the  optimum  spacing  is  one  of  the  objectives  of  this 
long-term  research  program. 

C.  CONDENSATION  RESEARCH  AT  NAVAL  POSTGRADUATE  SCHOOL 

The  research  effort  at  NPS  has  included  the  study  of  differing  fin  dimensions 
(i.e.  fin  height,  fin  width,  fin  spacing)  on  low-integral  finned  horizontal  tubes. 
Experimentation  has  included  the  use  of  three  difierent  test  fluids  (steam,  R-113, 
and  ethylene  glycol)  under  various  operational  conditions  using  a  number  of  different 
tube  diameters. 

Van  Petten  [Ref.  2]  provides  a  summary  of  research  efforts  at  NPS  through  the 
end  of  1988.  Van  Petten  and  subsequent  researchers  have  analyzed  small,  medium, 
and  large  diameter  finned  tubes  to  find  the  optimum  fin  spacing  for  maximum  heat- 
transfer  enhancement  of  the  fluids  mentioned  above.  However,  discrepancies  found 
by  Guttendorf  [Ref.  3]  in  the  data  processing  technique  (modified  Wilson  plot), 
which  resulted  in  different  values  of  heat-transfer  enhancement  (for  the  same  tube 
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under  the  same  operating  conditions  depending  on  whether  an  insert  was  or  was  not 
used),  have  raised  doubts  about  the  accuracy  of  the  inside  heat-transfer  correlation. 

Rouk  [Ref.  4]  investigated  the  use  of  an  optimization  technique  to  predict  the 
inside  heat-transfer  correlation.  When  the  optimization  effort  proved  unsuccessful, 
he  next  used  the  instrumented  smooth  tube  data  of  Georgiadis  [Ref.  5]  to  develop 
an  inside  heat  transfer  correlation,  but  could  not  find  a  correlation  with  sufficient 
accuracy  based  on  previous  data.  He  recommended  that  once  an  overhaul  on  the 
test  apparatus  was  complete,  an  increase  in  data  precision  would  allow  the 
development  of  an  accurate  inside  heat-transfer  correlation.  This  work  is  a  follow 
on  effort  to  develop  inside  heat-transfer  correlations  which  can  predict  the  value  of 
the  inside  heat-transfer  coefficient  with  good  accuracy  under  a  variety  of  flow 
conditions.  Once  an  inside  correlation  can  be  found,  the  object  of  this  effort  is  to 
reprocess  previous  data  and  see  if  the  discrepancies  reported  by  former  researchers 
for  finned  tubes  on  this  apparatus  can  be  rectified. 

D.  OBJECTIVES 

The  main  objectives  of  this  thesis  were  to: 

1.  Disassemble  and  meticulously  clean  the  apparatus  to  eliminate  any  existing 
contamination  with  a  view  to  eliminating  dropwise  condensation  problems 
experienced  in  the  past. 

2.  Carefully  reassemble  the  apparatus  using  new  gasket  material,  and  make 
modifications  to  improve  system  performance. 

3.  Recalibrate  all  system  instrumentation  to  ensure  the  greatest  achievable 
accuracy. 
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4.  Investigate  the  possibility  of  manufacturing  large,  medium,  and  small  diameter 
instrumented  smooth  tubes. 

5.  Use  the  new  instrumented  tubes  and  the  one  existing  medium  diameter 
smooth  instrumented  tube  (fabricated  by  Poole  [Ref.  6])  to  obtain  accurate 
inside  heat-transfer  correlations  for  a  number  of  insert  types  as  well  as  the 
no  insert  condition. 

6.  Evaluate  the  accuracy  of  the  currently  used  data  processing  technique 
(modified  Wilson  plot)  using  instrumented  tube  data. 

7.  Reprocess  previous  data  using  the  new  correlations  with  a  view  to  comparing 
current  and  past  smooth  tube  and  finned  tube  data  to  provide  continuity  with 
previous  studies. 


II.  LITERATURE  SURVEY 


A.  INTRODUCTION 

When  a  vapor  condenses  in  the  fUmwise  mode  on  a  smooth  horizontal  tube  it 
forms  a  thin  continuous  film  of  condensate  on  the  surface  of  the  tube.  The 
condensate  film  thickens  around  the  tube  due  to  gravity.  This  condensate  film 
provides  a  resistance  to  heat  transfer  which  may  be  lowered  through  the  use  of  fins. 
For  quite  some  time,  it  was  thought  to  be  impractical  to  use  finned  tubes  with  high 
surface  tension  fluids  such  as  water,  due  to  condensate  retention  and  flooding 
between  the  fins.  However,  a  number  of  studies  conducted  on  finned  tubes  using 
steam  have  shown  that  substantial  heat-transfer  enhancement  may  be  achieved. 

A  significant  amount  of  research  at  the  Naval  Postgraduate  School  and 
elsewhere  has  addressed  the  issue  of  optimum  fin  height,  thickness,  and  spacing 
required  for  maximum  heat  transfer.  Yau  et  al  [Ref.  7]  reported  that  "with  an 
increase  in  fin  density,  up  to  a  limit  (this  limit  is  not  yet  known  in  a  generalized 
manner),  the  heat-transfer  coefficient  increases  at  a  rate  faster  than  the  increase  in 
the  outside  area  due  to  the  presence  of  fins.  This  additional  enhancement  is  due  to 
the  thinning  effect  of  the  surface-tension  forces  on  the  condensate  film. 
Unfortunately,  surface-tension  forces  also  adversely  affect  heat  transfer  by  causing 
condensate  to  be  retained  between  fins"  [Ref.  8].  Katz  et  al  [Ref.  9]  also  found  that 
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on  finned  tubes  the  portion  of  the  surfeice  occupied  by  condensate  is  dependent  upon 
the  ratio  of  condensate  surface-tension  to  density  and  the  fin  geometry. 

Condensate  retention  and  the  behavior  of  the  condensate  film  on  the  tube 
surface  under  various  conditions  are  critical  parameters  in  the  heat  transfer  process 
on  horizontal  finned  tubes.  Several  models  have  been  developed  to  predict  this 
behavior  and  the  reader  is  referred  to  an  extensive  review  of  horizontal  finned  tube 
heat  transfer  by  Marto  [Ref.  10]  for  a  more  detailed  coverage  of  the  topic. 

B.  VAPOR  SIDE  CONSIDERATIONS 

The  filmwise  condensation  of  vapor  on  a  horizontal  tube  is  a  complex  two- 
phase  heat  transfer  process,  for  which  a  suitabfy  complex  model  would  be  required 
to  accurately  predict  heat  transfer  performance  under  aU  conditions. 

In  1916,  Nusselt  [Ref.  11]  set  forth  his  theoretical  work  on  the  study  of  laminar 
filmwise  condensation  of  a  "stationary"  vapor  on  a  vertical  or  inclined  plate  and  a 
horizontal  tube.  Nusselt’s  simplifying  assumptions  included  the  following  [Ref.  12]: 

1.  Pure  saturated  vapor 

2.  Negligible  vapor  velocity  (U,=0) 

3.  Heat  transfer  across  the  condensate  film  by  conduction  only 

4.  Laminar  condensate  flow  governed  only  by  gravitational  and  viscous  forces 

5.  Condensate  properties  constant 

6.  Isothermal  condensing  surface 

7.  Negligft}le  interface  temperature  drop 
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Nusselt’s  result  for  the  mean  heat-transfer  coefficient  for  a  horizontal  tube  was 


obtained: 


=0.728 


^P/P/-P>»A 


(2.1) 


or 


Afu=  0.655 


*/p/p/-p^gAy^ 


1/3 


(2.2) 


where: 

Nu  =  mean  Nusselt  number 

kf  =  thermal  conductivity  of  condensate  film  (W/m  •  k) 

Pi  =  condensate  film  density  (kg/m^) 

=  vapor  density  (kg/m^) 
g  =  gravitational  constant  (9.81  m/s^) 

h(g  =  specific  enthalpy  of  vaporization  (J/kg) 

Pi  =  dynamic  viscosity  of  condensate  film  (N  s/m^) 

Do  =  outside  tube  diameter  (m) 

ATf  =  average  temperature  difference  across  condensate  film  (K) 
q  =  heat  flux  based  on  outside  area  (Q/A^)  (W/m^) 


Many  workers  have  improved  on  Nusselt’s  theoretical  analysis  by  accounting 
for  some  of  the  terms  he  neglected  through  his  simplifying  assumptions.  However, 
equations  (2.1)  and  (2.2)  have  been  found  to  be  remarkably  accurate  over  a  wide 

range  of  conditions  for  a  stationary  vapor.  High  vapor  velocity  can  increase  film 

« 

condensation  heat  transfer  substantially.  This  enhancement,  which  refers  to  the 
amount  of  heat  transfer  above  or  below  the  value  predicted  by  the  Nusselt  analysis, 
is  due  to  the  effect  of  thinning  the  condensate  film.  However,  vapor  shear  is  the  one 
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assumption  which  if  applied  can  lead  to  significant  increases  in  the  heat-transfer 
coefficient. 

The  theoretical  result  of  Shekriladze  and  Gomelauri  (1966)  [Ref.  13],  who 
considered  interfacial  shear  stress  due  to  vapor  velocity,  is  shown  in  equation  (2.3). 


Nu 


=0.64(1 +(l  +  1.69f)‘^*^ 


(2.3) 


where: 

Nu  =  Nusselt  number  for  the  vapor  side 
Re24,  =  two  phase  Reynolds  number,  (py  U,  D^py) 


For  steam  condensation,  the  empiricaUy  derived  correlation  of  Fujii  et  al  [Ref. 
14]  is  shown  in  equation  (2.4). 


Nu 


=0.96FV* 


(2.4) 


The  Nusselt  expression  (equation  (2.1)  can  be  expressed  in  similar  form: 


Nu 

ReS 


=0J2»F^'* 


(2.5) 


Whereas  the  vapor  velocity,  U«,  cancels  out  in  the  Nusselt  expression 
(stationary  vapor  assumption),  the  Fujii  correlation  includes  the  vapor  velocity  effect 
Therefore  we  can  expect  equation  (2.4)  to  more  accurately  predict  steam-side  heat 
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transfer  coefficients  for  those  cases  where  vapor  velocity  begins  to  have  a  significant 
impact. 


For  further  review  of  basic  theoretical  studies  on  the  subject  of  laminar  film 
condensation  on  smooth  tubes  the  reader  is  referred  to  Rose  [Ref.  15]. 


C.  COOLANT  SIDE  CONSIDERATIONS 

For  a  turbulent  flow  regime  inside  a  pipe  (Re  >  10,000),  a  number  of  coolant- 
side  correlations  have  been  used;  many  of  these  have  taken  the  form: 


(2.6) 


where: 

Nu  =  mean  coolant  Nusselt  number  for  turbulent  flow 
C,  =  correlating  coefficient 
Re  =  coolant  Reynolds  number 
Pr  =  coolant  Prandtl  number 


The  most  common  correlations  with  the  same  form  as  equation  (2.6)  are  that 
of  Dittus  and  Boelter  (1930)  [Ref.  16]: 

Nu  =0.023  Re^*Pr^*  (2  '7) 


and  Colburn  (1933)  [Ref.  17]: 


Nu  =0.023  Re^*Pr^ 


(2.8) 
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Sieder  and  Tate  (1936)  [Ref.  18]  applied  a  correction  factor  to  equation  (2.8) 
to  account  for  cases  in  which  the  bulk  to  inner  wall  temperature  difference  is  large 
enough  to  cause  substantial  variations  in  coolant  viscosity  as  follows: 

Nu  =0.027  (210) 

where: 

=  coolant  viscosity  evaluated  at  mean  bulk  temperature  (N  s/m^) 

IJL^  =  coolant  viscosity  evaluated  at  mean  inner  tube  wall  temperature 
(Ns/m^) 

The  fluid  properties  in  equations  (2.7),  (2.8),  and  (2.9)  are  evaluated  at  mean 
coolant  bulk  temperature  T„, 

(211) 

where: 

T,  =  tube  coolant  inlet  temperature  (K) 

Tj  =  tube  coolant  outlet  temperature  (K) 

Equations  (2.7),  (2.8)  and  (2.9)  are  valid  for  Re>  10^  and  0.7  <  Pr  <  100,  and 
were  developed  for  long  smooth  pipes  with  no  inserts  [Ref.  12]. 

The  use  of  inserts  and  the  effect  of  bends  close  to  the  tube  entrance  region 
can  affect  the  values  of  both  the  leading  coefficient,  C,  ,  and  Reynolds  number 
exponent,  m,  in  equation  (2.6).  One  of  the  major  focuses  of  this  study  is  to 
determine  the  values  of  Q  and  m  for  the  no  insert,  wire  wrap  insert,  and  Heatex 
insert  cases. 
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Other  well-known  turbulent  pipe  flow  heat-transfer  correlations  (i.e.  Petukhov- 
Popov  [Ref.  19],  Sleicher-Rouse  [Ref.  20],  etc.)  and  the  results  of  an  ANL  (Argonne 
National  Laboratory)  [Ref.  21]  study  which  evaluated  several  such  correlations  for 
accuracy  are  reviewed  in  section  VI  C. 
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III.  APPARATUS  AND  SYSTEM  INSTRUMENTATION 

A.  SYSTEM  OVERVIEW 

The  apparatus  used  for  this  research  was  basically  the  same  as  was  used  by  Van 
Petten  [Ref.  2]  and  Guttendorf  [Ref.  3]  with  certain  modifications.  A  ^tem 
schematic  is  provided  in  Figure  1.  Steam  generated  from  distilled  water  in  the  .30 
m  diameter  pyrex  glass  boiler  using  ten  4  kW,  440  V  Watlow  immersion  heaters  was 
the  working  medium  for  this  set  of  experiments.  From  the  boiler  section  the  steam 
passed  up  through  a  reducing  section  and  a  2.13  m  straight  length  of  pyrex  glass 
piping,  (ID  of  0.15  m),  it  was  then  turned  through  180  degrees  using  two  90  degree 
pyrex  glass  elbows,  and  then  descended  down  a  1.52  m  straight  length  of  pyrex  glass 
piping.  The  steam  then  entered  the  stainless  steel  test  section  containing  the 
horizontally  mounted  condenser  tube  (see  Figures  1  and  2);  any  steam  not 
condensing  there  was  condensed  in  the  auxiliary  condenser  located  just  beneath  the 
test  section.  The  auxiliary  condenser  was  constructed  of  a  single  copper  coil  mounted 
to  a  stainless  steel  baseplate  enclosed  within  a  pyrex  glass  condenser  section. 
Coolant  flow  through  the  auxiliary  condenser  was  used  to  control  ^tem  pressure, 
and  all  condensed  liquid  was  returned  via  the  condenser  baseplate  drain  to  the  boiler 
section  by  gravity. 

Coolant  for  the  auxiliary  condenser  was  provided  via  a  throttled  water 
connection  with  associated  flowmeter.  Coolant  flow  through  the  single  horizontal 
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Figure  1.  Schematic  of  Single  Tube  Test  Apparatus 
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DRAIN  FILL /DRAIN  VALVE 


Figure  2.  Schematic  of  Test  Section  Insert 
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tube  was  provided  by  a  separate  system  consisting  of  a  sump  tank  with  two 
centrifugal  pumps  connected  in  series.  Coolant  flow  rate  was  measured  by  a 
carefully  calibrated  flowmeter.  By  varying  coolant  flow-rate  through  the  single 
horizontal  tube,  the  rate  of  steam  condensation  on  the  tube  (and  hence  heat-transfer 
coefficient)  could  be  varied. 

Non-condensible  gases  were  removed  using  the  vacuum  pump  system  shown  in 
Figure  3.  The  condensing  coil  for  this  purge  system,  located  in  the  sump  tank, 
served  to  condense  steam  carried  through  the  vacuum  line  during  the  purging 
process.  The  vacuum  line  took  its  suction  from  the  base  of  the  auxiliary  condenser, 
the  coolest  spot  in  the  apparatus  and  the  place  where  non-condensible  gases  (i.e.  air) 
were  most  likely  to  accumulate. 

B.  SYSTEM  INSTRUMENTATION 

The  power  to  the  440  V  heaters  was  controlled  through  a  panel  mounted 
potentiometer.  A  description  of  the  power  calculation  for  input  into  the  data 
acquisition  system  can  be  found  in  Poole  [Ref.  6]. 

System  pressure  was  monitored  in  three  ways: 

1.  A  Setra  model  204  pressure  transducer 

2.  A  Heise  solid  front  pressure  gauge  (visual  reading  only) 

3.  System  saturation  pressure  from  vapor  temperature  measurement 
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System  vapor  temperature  was  monitored  using  both  teflon,  and  metal  sheathed 
type-T  copper/constantan  thermocouples  located  juxtaposed  in  the  test  section;  this 
position  was  just  upstream  of  the  test  condenser  tube.  Condensate  temperature  was 
also  monitored  using  a  teflon  coated  type-T  copper/constantan  thermocouple  located 
on  the  condensate  return  line  between  the  auxiliary  condenser  and  boiler.  Coolant 
temperature  rise  in  the  condenser  tube  was  measured  using  four  methods: 

1.  Two  teflon  coated  type-T  copper/constantan  thermocouples 

2.  Two  metal  sheathed  type-T  copper/constantan  thermocouples 

3.  Two  Hewlett-Packard  2804A  quartz  crystal  thermometers 

4.  A  ten-junction  teflon  coated  type-T  copper/constantan  thermopile 

These  were  all  placed  at  the  inlet  to  and  exit  from  the  condenser  tube;  at  the 
outlet,  all  thermocouples  were  placed  just  downstream  of  a  coolant  mixing  chamber. 

Two  data  reduction  programs  were  used  to  collect  and  reduce  data  on  this 
apparatus;  "DRPINST',  and  "DRPKS".  The  instrumented  tube  constructed  by  Poole 
[Ref.  6]  was  used  to  determine  an  accurate  inside  heat-transfer  correlation  for  inserts 
used;  this  instrumented  tube  contained  six  wall  thermocouples.  For  the  instrumented 
tube  the  appropriate  calibration  equations  were  accessed  in  the  data  acquisition 
program  "DRPINST'.  For  non-instrumented  tubes,  the  data  reduction  program 
"DRPKS"  was  used.  Fluid  property  equations  used  in  the  data  reduction  programs 
are  given  in  Appendix  A  and  calibrations  were  conducted  for  all  system 
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instrumentation  (flowmeter,  thermocouples,  pressure  transducer,  etc.)  and  are 
included  in  Appendix  B. 

The  data  as  monitored  by  the  aforementioned  system  instrumentation  was 
processed  by  an  HP-3497A  data  acquisition  system  controlled  by  an  HP-9826A 
computer  provided  with  the  correct  data  acquisition  program.  The  raw  data  was 
processed  and  stored  on  computer  disks.  Program  channel  assignments  are  given  in 
Table  1. 

C.  SYSTEM  MODinCATIONS 

At  the  beginning  of  this  investigation  the  apparatus  was  entirely  disassembled 
to  facilitate  complete  overhaul  of  the  system.  Modifications  and  details  of  assembly 
were  as  follows: 


1 .  The  apparatus  was  taken  apart  piece  by  piece,  inspected,  meticulously  cleaned 
with  a  warm  solution  of  Sparkleen  biodegradable  soap  and  subjected  to  a 
complete  acetone  rinse  prior  to  reassembly. 

2.  A  new  pyrex  glass  riser  section  above  the  boiler  was  0.31  m  shorter  than  the 
previous  section  and  allowed  the  addition  of  a  new  aluminum  stand  on  which 
to  place  the  heater  baseplate.  This  new  stand  allowed  much  easier  access  to 
the  440  V  heater  wiring  plus  the  adjustable  legs  allowed  level  adjustment  of 
the  entire  apparatus,  to  ensure  proper  alignment  of  the  single  horizontal  tube. 

3.  The  two  pyrex  glass  elbows  were  replaced. 

4.  Every  gasket  in  the  system  was  replaced  (using  Buna-N  rubber)  and,  using  a 
standard  star  torque  pattern,  all  flanged  joints  were  tightened  to  a  final 
torqued  of  30  inch-pounds  (manufacturer  recommended  maximum  torque 
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Table  1.  DATA  ACQUISITION  SYSTEM  CHANNEL  ASSIGNMENT 


CHANNEL 

QUANTITY 

APPLICABLE  TYPE 

MEASURED 

DESIGNATION 

40 

Vapor  Temperature 

T-57 

(See  Appendix  B.l) 

41 

Vapor  Temperature 

T-56 

(See  Appendix  B.l) 

42 

Room  Temperature 

T-58 

(See  Appendix  B.l) 

43 

Tube  Coolant  in  1 

T-58 

44 

Tube  Coolant  out  1 

T-58 

45 

Tube  Coolant  in  2 

T-55 

(See  Appendix  B.l) 

46 

Tube  Coolant  out  2 

T-55 

47 

Condensate  Return 

T-58 

48 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

49 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

50 

Instrumented  Tube 

T-57  or  T-55 
as  appUcable 

51 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

52 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

53 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

60 

10  -  Junction 
Thermopile 

T-61 

(See  Appendix  B.l) 

61 

« 

Voltage 

NA 

62 

Current 

NA 

64 

Pressure  Transducer 

(See  Appendix  B.4) 
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specification  was  60  inch-pounds).  The  bolts  holding  the  flanges  together 
should  be  checked  for  tightness  on  a  regular  basis,  since  the  thermal  cycling 
of  the  apparatus  has  been  shown  to  cause  loosening  of  nut  and  bolt 
assemblies. 

5.  The  previous  vacuum  pump  was  replaced  with  a  Gast  model  2567-V108 
vacuum  pump  which  could  draw  vacuum  to  130  mmHg  with  an  installed  check 
valve  to  prevent  pump  back  spin  when  the  pump  was  stopped.  The  new 
pump  could  draw  vacuum  much  more  rapidly,  but  could  not  operate  at  less 
than  130  mmHg.  Remaining  non-condensible  gases  were  removed  by  flushing 
the  system  with  steam.  The  steam  flushing  procedure  for  removal  of  non¬ 
condensible  gases  is  given  in  the  operating  procedures  section.  Once  the 
system  is  completely  filled  with  steam,  operating  pressures  well  below  130 
mmHg  could  be  acUeved  utilizing  the  amdliary  condenser. 

6.  The  double  coil  auxiliary  condenser  was  replaced  with  a  single  coil  used 
originally  by  Van  Petten  [Ref.  2].  The  single  coil  was  not  coated  with  the 
special  oxide  coating  used  by  Guttendorf  [Ref.  3].  It  was  felt  to  be  superior 
to  the  double  coil  in  that  the  baseplate  welds  were  of  much  higher  quality  and 
were  preferred  on  the  basis  of  vacuum  tightness. 

7.  The  aluminum  side  plates  attached  to  the  pyrex  glass  auxiliary  condenser 
housing  were  replaced  with  new  stainless  steel  side  plates  with  penetrations 
for  pressure  bleed,  vacuum  line,  and  a  pressure  transducer.  These  three 
penetrations  were  fitted  with  screw  threaded  stainless  steel  connectors.  The 
stainless  steel  connectors  were  heli-arc  welded  in  place.  Prior  to  this 
modification  (completed  24  January  1992)  a  leak  test  conducted  from  21 
December  1991  to  2  January  1992  revealed  a  mean  vacuum  leak  rate  of  ~  3.4 
mmHg  per  day  (see  Figure  C.l).  A  subsequent  leak  test  conducted  fi-om  6 
February  1992  to  19  February  1992  showed  an  improvement  in  the  mean  leak 
rate  to  —1.7  mmHg  per  day  (see  Figure  C.2). 

8.  System  instrument  modifications  included  the  addition  of  the  Setra  pressure 
transducer  and  the  Heise  pressure  gauge,  and  the  removal  of  the  mercury 
manometer.  All  system  instrumentation  was  recalibrated  and  the  results 
incorporated  into  the  data  reduction  programs. 

9.  Finally  the  apparatus  was  lagged  with  Halstead  insulating  foam  to  reduce  heat 
loss  as  much  as  possible.  The  test  section,  which  was  left  uncovered 
previously,  was  also  lagged. 


22 


IV  EQUIPMENT  OPERATION  AND  EXPERIMENTAL  PROCEDURE 


A.  SYSTEM  STARTUP  AND  SHUTDOWN  PROCEDURES 

Startup  of  the  system  is  accomplished  in  the  following  manner: 

1.  Ensure  distilled  water  level  in  the  boiler  is  4  to  6  inches  above  the  top  of  the 
heating  elements.  The  boiler  is  filled  by  gravity  drain  via  a  hose  connection 
from  the  distilled  water  tank  to  the  boiler  fill  valve.  Ensure  the  vent  valve  on 
the  side  of  the  auxiliary  condenser  is  open  when  filhng  or  draining  the  boiler. 
The  boiler  may  be  drained  by  removing  the  hose  connection  and  opening  the 
fill  valve,  which  allows  drainage  into  the  trench  directly  beneath  the  boiler. 

2.  Once  the  boiler  is  filled  to  the  appropriate  level,  shut  the  boiler  fill  valve  and 
the  distilled  water  tank  valve. 

3.  Shut  the  system  vent  valve. 

4.  Turn  on  the  data  acquisition  system,  computer  and  printer.  Load  the 
appropriate  program  (DRPINST  or  DRPKS)  and  check  for  proper  operation. 
TTien  check  all  thermocouple  outputs,  by  stepping  through  the  appropriate 
data  acquisition  system  channels,  to  verify  that  all  are  registering  ambient 
temperature. 

5.  Open  the  fill  valve  to  the  coolant  water  sump  tank  to  a  level  such  that  the 
tank  overflow  drain  box  does  not  overflow  (the  valve  is  located  between  the 
boiler  control  panel  and  heat  pipe  apparatus). 

6.  Turn  on  the  cooUng  water  supply  pumps  and  adjust  the  tube  flow  rate  from 
20%  to  60%  of  the  rotameter  setting  and  check  for  leaks.  Reset  flow  rate  to 
desired  level. 

7.  Open  valves  from  tap  water  system  to  auxiliary  condenser  and  adjust  coolant 
flow  rate  to  at  least  30%  and  check  for  leaks.  Reset  flow  rate  to  at  least 
10%. 

8.  Energize  heaters  and  adjust  voltage  to  approximately  50  volts  (40  volts  if  the 
system  is  already  at  vacuum  below  100  mmHg  to  limit  the  vibrational  shock 
to  the  system  from  oversized  vapor  bubble  formation).  To  energize  the 
heaters  there  are  three  switches  which  must  be  placed  in  the  on  position.  The 
first  is  located  in  power  panel  p5  located  in  the  main  hallway  adjacent  to  the 
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lab  and  is  labeled  switch  3  /  heater  controUer  room  106.  The  second  is  the 
heater  load  bank  circuit  breaker  located  on  the  side  of  the  boiler  control 
panel.  The  third  is  the  condensing  rig  boiler  power  switch  on  the  front  of  the 
boiler  control  panel.  Increase  the  voltage  gradual^  in  10  volt  increments  to 
the  desired  level. 

9.  Turn  on  the  vacuum  pump  and  open  the  vacuum  line  valve.  Allow  the 
vacuum  pump  to  run  until  system  pressure  is  below  3  psi,  then  shut  the 
vacuum  line  valve  just  prior  to  turning  ofr  the  vacuum  pump. 

10.  As  system  warmup  continues  and  pressure  increases  to  above  4  psi  energize 
the  vacuum  pump  as  necessary  to  flush  the  non-condensible  gases  out  of  the 
system  through  the  vacuum  line  by  forcing  the  gases  out  with  steam.  To 
ensure  that  the  non-condensibles  gather  at  the  base  of  the  auxiliary 
condenser,  where  the  vacuum  line  suction  is  located,  ensure  that  the 
horizontal  tube  is  not  supplied  with  coolant  flow,  and  adjust  coolant  flow 
through  the  auxiliary  condenser  as  necessary  to  ensure  steam  is  filling  the 
entire  system.  The  auxiliary  condenser  may  be  touched  lightly  by  hand,  along 
its  entire  length,  to  ensure  the  system  is  completely  filled  with  steam;  any  cool 
spots  indicate  the  presence  of  non-condensible  gases  which  means  that  the 
flushing  process  is  not  complete.  The  flushing  process  takes  15  to  30  minutes 
to  accomplish,  and  should  be  repeated  periodically  for  long  periods  of 
operation. 

11.  At  the  conclusion  of  the  flushing  process,  shut  the  vacuum  line  valve  and 
secure  the  vacuum  pump. 

12.  In  order  to  ensure  that  filmwise  condensation  occurs  on  the  tube,  coolant  flow 
through  the  tube  must  be  initiated  as  follows: 


a.  Allow  the  apparatus  vapor  temperature  measurement  (channel  40)  to 
reach  at  least  3800  microvolts. 

b.  Cut  in  the  auxiliary  coolant  flow  (50%  or  60%  level)  to  cool  the  vapor 
temperature  to  roughly  3200  microvolts. 

c.  Secure  coolant  flow  through  the  auxiliary  condenser,  and  allow  the 
vapor  temperature  level  to  climb  to  about  3700-3800  microvolts,. which 
allows  a  steam  blanket  to  cover  the  tube. 

d.  Initiate  coolant  flow  through  the  single  horizontal  tube  at  the  80%  level. 
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e.  Cut  in  coolant  flow  to  the  auxiliary  condenser  to  control  pressure,  and 
observe  the  condensation  process  to  ensure  that  a  condensate  film  has 
formed  on  the  tube. 

13.  Run  software  program  DRPINST  for  an  instrumented  tube  (DRPKS  for 
uninstrumented  tubes)  by  pressing  "run"  on  the  keyboard. 

To  take  data  for  an  instrumented  tube,  the  questions  for  DRPINST  can  be 

answered  as  follows: 

•  Select  fluid  ...  Enter  0  for  steam 

•  Select  option  ...  Enter  1  to  take  new  data 

•  Enter  month,  date  and  time  ...  Press  enter 

•  Enter  input  mode  ...  Enter  1  for  new  data 

•  Give  a  name  for  the  raw  data  file  ...  Enter  name 

•  Enter  geometry  code  ...  Enter  1  for  finned,  0  for  plain 

•  Select  insert  type  ...  Enter  0  for  none,  1  for  twisted  tape,  2  for  wire  wrap,  3  for 

Heatex 

•  No.  of  thermocouples  in  wall?  ...  Enter  4,  5,  or  6  depending  on  the  tube 

•  Select  tube  diameter  type  ...  Enter  2  for  medium 

•  Enter  pressure  condition  ...  Enter  0  for  vacuum,  1  for  atmospheric 

•  Give  a  name  for  the  wall  temperature  file  ...  Enter  name 

•  Select  input ...  Enter  1  for  short,  2  for  long,  or  3  for  raw  data 

•  Like  to  check  NG  (non-condensible  gas)  concentration ...  Enter  1  for  yes,  2  for 

no;  you  must  answer  yes  for  the  first  data  point 

•  Enter  flowmeter  reading  (%)  ...  Enter  2  digit  number  (i.e.  20  or  58  etc.) 

•  Connect  voltage  line  ...  Flip  the  voltage  line  toggle  switch,  located  on  the  power 

control  panel,  to  the  on  position  and  press  enter 
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•  Disconnect  voltage  line ...  Flip  the  voltage  line  toggle  switch  off  and  press  enter 

•  Enter  pressure  gauge  reading  (Pga) ...  Enter  reading  off  gauge  in  psi 

•  Select  measurement ...  Enter  0  for  teflon,  1  for  metal  sheath,  2  for  quartz, 

3  for  thermopile 

•  Change  TC(X)L  rise?  ...  Enter  1  for  yes,  2  for  no 

•  OK  to  store  this  data  set?  ...  Enter  1  for  yes,  0  for  no 

•  Will  there  be  another  run?  ...  Enter  1  for  yes,  0  for  no;  starts  at  check  NG 

concentration  for  following  runs. 

To  take  data  for  an  uninstrumented  tube  the  questions  for  DRPKS  can  be 
answered  as  follows: 

•  Select  fluid  ...  Enter  0  for  steam 

•  Select  option  ...  Enter  0  to  take  new  data 

•  Enter  month,  date  and  time  ...  Press  enter 

•  Enter  disk  number  ...  Enter  number 

•  Enter  input  mode  ...  Enter  0  for  new  data 

•  Select  Ci ...  Enter  0  to  find  a  Ci  value,  2  to  use  a  Ci  value  stored  in  the 

program. 

•  Give  a  name  for  the  raw  data  file  ...  Enter  name 

•  Enter  geometry  code  ...  Enter  1  for  finned,  0  for  plain 

•  Enter  insert  type  ...  Enter  0  for  none,  1  for  twisted  type,  2  for  wire  wrap,  3  for 

Heatex 

•  Select  tube  type  ...  Enter  0  for  thick  waU  (only  thick  wall  tubes  were  tested) 

•  Select  material  code  ...  Enter  0  for  copper  (only  copper  tubes  were  tested) 
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•  Select  tube  diameter  type  ...  Enter  1  for  medium  (no  small  or  large  diameter 

tubes  were  tested) 

•  Enter  pressure  condition  ...  Enter  0  for  vacuum,  1  for  atmospheric 

•  Want  to  create  a  file  for  NR  vs  F?  ...  Enter  1  for  yes,  0  for  no 

•  Give  a  name  for  plot  data  file  ...  Enter  name;  easiest  to  use  the  raw  data  file 

name  preceded  by  a  P 

•  Select  output ...  Enter  0  for  short,  1  for  long,  2  for  raw  data 

•  Like  to  check  NG  concentration  ...  Enter  1  for  yes,  2  for  no;  you  must  answer 

yes  for  the  first  data  point 

•  Enter  flowmeter  reading  (%)  ...  Enter  2  digit  number  (i.e.  20  or  60  etc.) 

•  Connect  voltage  line  ...  Flip  the  voltage  line  toggle  switch  on  and  press  enter 

•  Disconnect  voltage  line  ...  Flip  the  voltage  line  toggle  switch  off  and  press  enter 

•  Enter  pressure  gauge  reading  (Pga) ...  Enter  reading  off  gauge  in  psi 

•  Select  measurement  ...  Enter  0  for  teflon,  1  for  metal  sheath,  2  for  quartz, 

3  for  thermopile 

•  Change  TCOOL  rise  ...  Enter  1  for  yes,  2  for  no 

•  OK  to  store  this  data  set ...  Enter  1  for  yes,  0  for  no 

•  Will  there  be  another  run  ...  Enter  1  for  yes,  0  for  no;  starts  at  check  NG 

concentration  for  following  runs. 

14.  Only  answer  the  program  questions  up  to  "Enter  flowmeter  readings". 
Monitor  system  temperature  using  the  vapor  thermocouple  voltage  reading 
(the  program  automatically  resets  to  channel  40)  closely  until  system  warmup 
is  complete. 

15.  Monitor  system  temperature  and  pressure  carefuUy  to  prevent  a  system  over 
pressure  during  warmup  (especially  at  atmospheric  conditions). 
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16.  If  conducting  a  vacuum  run,  gradually  adjust  voltage  to  90  volts  (usually  in  10 
volt  increments).  Obtain  the  desired  q>erating  condition  by  manually 
controlling  coolant  flow  through  the  auxiliary  condenser  until  channel  40 
reads  1970  ±  20  microvolts  (~48°C).  Vapor  velocity  ~  2  m/s. 

17.  If  conducting  an  atmospheric  run,  gradually  adjust  voltage  to  175  volts  from 
the  90  volt  level  in  10-20  volt  increments.  Again  the  desired  operating 
condition  is  obtained  by  manually  controlling  coolant  flow  through  the 
auxiliary  condenser  until  channel  40  reads  4280  ±  20  microvolts  (~100°C). 
Vapor  velocity  ~  1  m/s. 

18.  Monitor  the  condensation  process  using  the  glass  viewing  window  periodically 
to  ensure  that  filmwise  condensation  is  maintained.  To  clear  the  viewing 
window  of  fog  and  moisture  increase  coolant  flow  through  the  auxiliary 
condenser  briefly  to  50%  or  60%,  then  reset  to  desired  flow  rate. 

19.  When  taking  readings  be  sure  to  check  the  flowmeter  setting  prior  to  entering 
it  into  the  computer  (it  has  a  tendency  to  fluctuate  slightly). 

20.  If  conducting  vacuum  and  atmospheric  runs  on  the  same  day  always  conduct 
the  vacuum  run  first.  If  the  atmospheric  run  is  done  first  it  takes  too  long  for 
the  system  to  cool  down  to  vacuum  operating  temperatures. 

The  system  is  secured  in  the  following  manner: 

1.  Secure  power  to  the  heating  elements. 

2.  Secure  coolant  flow  through  the  tube,  through  the  auxiliary  condenser,  and 
to  the  sump  tank. 

3.  If  desired  to  maintain  the  system  at  vacuum  conditions  until  the  next  run  the 
shutdown  is  complete.  Continued  cooling  water  circulation  may  be  used  to 
assist  in  cooling  down  the  system. 

4.  To  bring  the  system  back  to  atmospheric  conditions  slowly  open  the  vent 
valve. 

5.  The  data  acquisition  system  may  be  turned  off  whenever  it  is  not  necessary 
to  monitor  system  parameters. 

6.  Periodically  change  distilled  water  in  the  boiler. 
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7.  If  an  emergency  should  arise  such  as  abrupt  overpressurization  or  breakage, 
immediately  secure  power  to  the  heaters  and  open  the  vent  valve,  then  let  the 
system  cool  down  before  checking  the  apparatus  for  damage. 

B.  EXPERIMENTAL  PROCEDURES  AND  OBSERVATIONS 

Water  is  a  poor  wetting  medium  and  therefore  great  care  was  taken  to  ensure 

that  uniform  filmwise  condensation  was  the  only  condensation  mode  occurring  during 

a  data  run.  Even  though  the  apparatus  was  meticulously  cleaned  (as  mentioned 

previously),  a  continuing  problem  with  dropwise  condensation  manifested  itself. 

Subsequent  to  steam  cleaning  the  system  with  a  Sparkleen  soap  solution,  by 

operating  the  system  with  a  soapy  solution  in  the  boiler  (the  solution  bubbled 

through  the  entire  apparatus),  dropwise  condensation  was  observed  on  the  installed 

instrumented  tube.  After  taking  some  data  when  in  the  dropwise  condition,  the  tube 

was  removed  and  rigorously  cleaned  using  a  warm  Alconox  soap  solution  with  a 

scrub  brush.  However,  after  observing  the  filmwise  mode  initially,  the  condensation 

mechanism  soon  transitioned  to  mixed  mode  and  then  back  to  the  dropwise  mode. 

Since  only  a  filmwise  condition  over  several  hours  would  suffice,  the  tube 

chemical  treatment  procedure  used  by  Guttendorf  [Ref.  3]  and  several  other 

researchers  was  used  to  produce  filmwise  condensation.  The  tube  was  chemically 

treated  prior  to  installation  as  follows: 

1.  Clean  the  internal  and  especially  the  external  surfaces  of  the  tube  using  a  soft 
brush  and  mild  soap  (using  the  Alconox  detergent  in  warm  water),  rinse  with 
acetone  then  rinse  thoroughly  with  distilled  water.  Repeat  the  cleaning 
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procedure  until  the  distilled  water  rinse  perfect^  wets  the  tube  surface;  any 
breaks  in  the  wetting  film  at  this  point  are  likely  to  result  in  dropwise 
condensation  spots  once  the  tube  is  instaUed  in  the  apparatus. 

2.  Place  the  tube  in  a  steam  bath. 

3.  Mix  equal  amounts  of  ethyl  alcohol  and  a  50%  by  weight  solution  of  sodium 
hydroxide.  Keep  the  solution  warm  so  that  a  watery  consistency  is 
maintained. 

4.  Apply  the  solution  to  the  tube  with  a  small  paint  brush,  retaining  the  tube  in 
the  steam  bath.  If  the  tube  has  not  been  treated  previously,  apply  a  coating 
of  the  solution  every  10  minutes  for  an  hour.  If  the  tube  has  been  previously 
treated,  apply  a  coating  every  5  minutes  for  a  period  of  20  minutes. 

5.  Remove  the  tube  from  the  steam  bath  and  thoroughly  rinse  the  tube  with 
distilled  water  to  remove  any  excess  solution.  Install  the  tube  in  the  test 
section  immediately,  being  careful  not  to  touch  the  tube  surface.  Oil  or  dirt 
from  any  source  may  contaminate  the  tube  surface  and  result  in  mixed  mode 
or  dropwise  condensation. 


The  oxide  layer  which  forms  ou  the  tube  is  very  thin,  and  has  negligible  thermal 
resistance  and  high  wetting  characteristics. 

Once  the  tube  was  installed  in  the  apparatus  (with  the  desired  insert  in  place), 
the  system  startup  procedure  outlined  in  section  IV  A  was  followed  to  take  data  at 
desired  conditions. 

At  vacuum  conditions,  when  single  tube  coolant  flow  was  initiated  with  vapor 
velocity  at  ~  2  m/s,  the  condensation  on  the  tube  did  not  develop  as  a  perfect  film 
but  instead  left  patches  where  the  film  was  broken.  These  patches,  or  streaks 
seemed  to  occur  at  regular  intervals,  and  it  was  postulated  that  they  were  due  to 
vortex  shedding  of  vapor  around  the  tube.  Therefore,  the  procedure  in  the  startup 
section  IV  A,  step  12,  was  used  to  promote  the  development  of  uniform  filmwise 
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condensation  by  inducing  a  stationary  vapor  condition  around  the  tube.  This  allowed 
a  steam  blanket  to  form  around  the  tube  prior  to  coolant  flow  initiation.  After  flow 
initiation  the  appearance  of  the  condensate  film  on  the  tube  surface  was  continuous 
with  no  breaks.  Momentary  film  instabilities  were  observed  at  vacuum  conditions  at 
pressures  below  ~  20  kPa  after  a  continuous  film  was  established  at  higher  pressures. 
These  instabilities  seemed  to  interrupt  the  film  sheet  only  for  an  instant  and  then 
disappear,  and  may  have  been  caused  by  vortex  shedding  of  the  vapor  around  the 
tube  as  already  mentioned.  A  possible  mechanism  to  explain  these  instabilities  is 
that  the  higher  vapor  velocity  at  vacuum  conditions  momentarily  thins  the  condensate 
film  via  vortex  shedding,  yet  this  thinning  effect  is  overcome  by  surface  tension  forces 
in  the  film  sheet  which  tend  to  restore  the  continuous  film.  These  instabilities  could 
only  be  seen  for  an  instant  and  then  would  vanish,  being  very  transitory  in  nature. 
There  also  did  not  appear  to  be  any  pattern  whatsoever  to  the  instability  formation. 
The  instabilities  were  not  observed  at  pressures  above  ~  20  kPa. 

The  data  taking  regimen  for  each  data  set  involved  starting,  then  verifying  the 
existence  of  a  filmwise  condensation  condition,  then  taking  data  at  flow  rates  (in  %) 
of  80,  70,  60,  50,  40,  30,  and  20  then  back  to  80  and  50  to  check  for  repeatability 
within  the  data  set.  Two  data  points  were  taken  at  each  of  the  first  seven  data  points 
and  one  each  for  the  last  two,  which  gave  a  total  of  16  data  points.  It  was  usually 
quite  clear  from  the  two  comparison  points,  and  fi’om  data  taken  previously  under 
similar  conditions,  whether  the  data  set  should  be  rejected  or  accepted.  After  tube 
installation,  the  appearance  of  one  or  more  small  patches  (breaks  in  the  film)  after 
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~  7-10  hours  of  operation  signaled  the  beginning  of  tube  contamination  which  got 
worse  with  time.  The  tube  would  then  be  removed  and  cleaned. 

C.  TUBES  TESTED 

The  data  taken  during  this  study  involved  extensive  use  of  the  instrumented 
smooth  tube  (SOI)  fabricated  by  Poole  [Ref.  6],  (six  wall  thermocouples  spaced  60 
degrees  apart  placed  at  midwall  and  midlength).  Due  to  excessive  thermocouple 
wear,  only  5  thermocouples  in  this  instrumented  tube  functioned  properly.  The  tube 
was  positioned  in  the  apparatus  to  make  optimum  use  of  functioning  thermocouples. 
The  preferred  arrangement  placed  4  thermocouples  at  10®,  190°,  250°,  and  310°  from 
the  top  dead  center  position  of  the  tube.  This  arrangement  provided  readings  from 
the  top  and  bottom  of  the  tube  and  two  intermediate  points,  giving  the  most  accurate 
mean  tube  wall  temperature  and  best  temperature  profile  readings.  At  the 
conclusion  of  this  study,  the  manufacture  of  the  new  instrumented  tubes  was  not  yet 
complete;  progress  to  date  (March  1992)  is  recorded  in  Appendix  F. 

Data  was  also  taken  on  a  uninstrumented  smooth  tube  (S02)  and  four  finned 
tubes  (803,  804,  80S,  806)  with  fin  spacings  of  0.5  mm,  1.0  mm,  and  1.5  mm,  and  2.0 
mm.  All  tubes  tested  were  classified  as  medium  tubes,  with  an  outside  diameter  of 
19.05  mm,  and  an  inside  diameter  of  12.70  mm.  The  four  finned  tubes  aU  had  the 
same  fin  height  of  1.0  mm,  and  fin  width  of  1.0  mm.  Data  runs  were  taken  either 
with  no  insert,  with  the  wire  wrap  or  the  Heatex  insert  installed.By  more  efficiently 
mixing  the  coolant  an  insert  significantly  increases  the  inside  heat-transfer  coefficient 
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The  Heatex  insert  consists  of  a  central  wire  core  onto  which  are  wound  a  series 
of  wire  loops,  each  inclined  at  a  common  angle  to  the  core.  The  loops  come  into 
direct  contact  with  the  tube  wall,  and  each  loop  provides  a  significant  amount  of 
coolant  mixing  as  the  coolant  flows  through  the  loop  mesh.  [Ref.  12]. 

The  wire  wrap  insert  was  a  copper  wire  spirally  wrapped  around  a  central 
stainless  steel  rod  with  a  uniform  pitch.  This  insert  induced  a  swirling  coolant 
motion,  which  enhanced  turbulent  mixing  within  the  tube.  This  particular  wire  wrap 
insert  was  the  same  insert  used  by  Guttendorf  [Ref.  3],  Coumbes  [Ref.  22],  and  Van 
Petten  [Ref.  2]  to  collect  data  on  the  medium  family  of  tubes. 

A  summary  of  data  runs  is  given  in  Table  2. 


Table  2.  DATA  RUN  SUMMARY 


Tube 

Filename 

Pressure 

(kPa) 

Vapor 

Velocity 

(m/s) 

Insert 

SOI 

DSOIMVSH2 

11 

2.17 

Heatex 

DSOIMASH4 

101 

1.06 

Heatex 

DSOIMASH5 

101 

1.06 

Heatex 

FSOIMASH3 

11 

2.11 

Heatex 

FIMAVSHl 

28 

3.53 

Heatex 

FIMAVSH2 

41 

2.45 

Heatex 

FIMAVSH3 

69 

1.53 

Heatex 

FIMAVSH4 

101 

1.08 

Heatex 

FIMAVSH5 

28 

1.01 

Heatex 

FIMAVSH6 

41 

1.04 

Heatex 

FIMAVSH7 

69 

1.02 

Heatex 

FIMASW3 

101 

1.07 

wire  wrap 

FIMASW4 

101 

1.07 

wire  wrap 
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Tube 


Fllenauie 


Insert 


FIMASW5 

COMPAIWl 

COMPAIW2 

FIMASN4 

FIMASN5 

FIMASN6 

COMPAINl 

COMPAIN2 

502  FNMAVSHl 
FNMAVSH2 
FNMAVSH3 
FNMAVSH4 
FNMAVSH5 
FNMAVSH6 
FNMAVSH7 
FNMAVSH8 

503  FIMAF051 
FIMAF052 

504  FIMAFlOl 
FIMAF102 

805  FIMAF151 

FIMAF152 


Pressure 

(kPa) 

Vapor 

Velocity 

(m/s) 

101 

1.06 

101 

1.07 

101 

1.07 

101 

1.08 

101 

1.07 

101 

1.08 

101 

1.08 

101 

1.08 

15 

6.24 

41 

2.42 

69 

1.51 

102 

1.03 

28 

1.02 

41 

1.03 

69 

1.02 

101 

1.07 

101 

1.06 

101 

1.06 

101 

1.07 

101 

1.06 

101 

1.07 

101 

1.07 

101 

1.07 

101 

1.07 

wire  wrap 
wire  wrap 
wire  wrap 
none 
none 
none 
none 
none 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 

Heatex 


FIMAF201 

FIMAF202 
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Tube  descriptions: 

501  instrumented  smooth  tube  fabricated  by  Poole 

502  uninstrumented  smooth  tube 

503  0.5  mm  finned  tube 

504  1.0  mm  finned  tube 

505  1.5  mm  finned  tube  (this  tube  designated  as  F096  by  Guttendorf;  also  a 
similar  1.5  mm  finned  tube  designated  as  F006  by  Guttendorf  and  Van 
Petten  was  not  testedl;  the  1.5  mm  finned  tube  results  were  not  included 
in  subsequent  finnea  tube  analysis  since  the  reason  for  differences  in 
experimental  results  for  these  two  tubes  has  not  yet  been  resolved. 

506  2.0  mm  finned  tube 
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V.  THEORETICAL  BACKGROUND  AND  DATA  REDUCTION  PROCEDURES 


A.  THEORETICAL  BACKGROUND 

The  overall  or  total  resistance  to  heat  transfer  from  vapor  to  coolant  consists 
of  the  sum  of  the  vapor  side  resistance  (R«),  the  tube  wall  resistance  (R„),  and  the 
coolant  side  resistance  (R,);  this  neglects  any  fouling  resistance  since  clean  tubes  are 
always  used. 


(5,1) 

The  vapor  and  coolant  side  resistances  are  convective  in  nature  and  may  be 
expressed  by  the  reciprocal  of  their  respective  heat-transfer  coefficient  and  surface 
area  product. 


(5.2) 


(5.3) 


where: 

R„  =  outside  vapor  side  resistance  to  heat  transfer  (K/W) 
ho  =  outside  heat-transfer  coefficient  (W/m^  K) 

Ao  =  effective  outside  surface  area  (m^) 

R,  =  inside  coolant  side  resistance  to  heat  transfer  (K/W) 
h;  =  inside  heat-transfer  coefficient  ^/m^*K) 

Af  =  effective  inside  surface  area  (m^) 

The  tube  wall  resistance  is  conductive  in  nature  and  is  represented  by  the  radial 
conduction  equation. 
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(5.4) 


l%Lk^ 


where: 

=  tube  wall  resistance  (K/W) 

Dr  =  outside  or  root  (if  finned)  diameter  (m) 

D,  =  inside  tube  diameter  (m) 

L  =  active  condensinjg  len^h  (133  mm) 

k,„  =  thermal  conductivity  of  tube  wall  (W/m-K) 

The  effective  outside  area  of  the  tube  is  calculated  using  the  following 
expression: 


=  (5.5) 

where: 

D,  =  outside  or  root  (if  finned)  diameter  (m) 

L  =  active  condensing  length  (133  mm) 

The  effective  inside  area  includes  the  inside  surface  area  involving  the  active 
condensing  length  and  the  inside  surface  area  of  the  insulated  inlet  and  outlet 
portions  of  the  tube.  These  portions  of  the  tube  act  as  fins  and  remove  heat  via 
axial  conduction.  The  extended  fin  assumption  with  associated  fin  efficiencies  was 
used  to  account  for  these  end  losses. 


A.  =  nD^(L  +  ^\^L^+r\2L2) 


(5.6) 


where: 

D,  =  inside  diameter  of  tube  (m) 

L  =  active  tube  condensing  length  (m) 

77i  =  fin  efficient^  of  inlet  portion  of  tube 

L]  =  length  of  inlet  portion  of  tube  (m) 

172  =  fin  efficiency  of  outlet  portion  of  tube 

Lj  =  length  of  outlet  portion  of  tube  (m) 

The  overall  thermal  resistance  to  heat  transfer  may  be  expressed  by: 
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(5.7) 


1 


Substituting  equations  (5.2)  and  (5.3)  into  equation  (5.7)  yields: 


1 

o  o 


K 


h,A, 


(5.8) 


where: 

U„  =  overall  heat-transfer  coefficient  (W/m^  K) 

Ao  =  effective  outside  surface  area  (nr) 
ho  =  outside  heat-transfer  coefficient  (W/m^  K) 

=  Tube  wall  resistance  (K/W) 
h,  =  inside  heat-transfer  coefficient  (WAn^  K) 

The  single  tube  condenser  apparatus  uses  the  log  mean  temperature  difference 
(LMTD)  analysis  for  calculation  of  the  heat  transfer  between  the  hot  vapor  and  cold 
coolant. 


Q^UJ,(lMTD)  (5.9) 

where: 

Q  =  heat  transfer  rate  to  the  cooling  water  (W) 

Uo  =  overall  heat-transfer  coefficient  (W/m^  K) 

=  effective  outside  surfiice  area  (m^) 

LMTD  =  log  mean  temperature  difference  between  vapor  and  coolant  (K) 
The  log  mean  temperature  difference  (LMTD)  is  given  by: 


LW7D  = 


di-r,) 


(5.10) 
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where; 


Ti  =  coolant  inlet  temperature  (K) 

T2  =  coolant  outlet  temperature  (K) 
Ts„  =  vapor  saturation  temperature  (K) 


In  this  and  previous  studies  at  NPS,  the  quartz  thermometer  ou^ut  for  T,  and 
Tj  were  used  in  the  calculations  for  the  coolant  temperature  rise,  and  the  saturation 
temperature,  T„„  was  measured  using  the  vapor  thermocouple  (channel  40). 

The  total  heat  transfer  across  the  tube  is  experimentally  determined  by 
measuring  the  mass  flow  rate  of  fluid  through  the  tube  and  its  accompanying 
temperature  rise. 

<?=mc/rj-r,)  (5.11) 

where: 

Q  =  heat  transfer  rate  (W) 
m  =  mass  flow  rate  of  coolant  (kg/s) 

Cp  =  specific  heat  of  coolant  at  constant  pressure  (J/kg-K) 

Tj  =  coolant  inlet  temperature  (K) 

T2  =  coolant  outlet  temperature  (K) 

Equation  (5.11)  may  be  used  directly  from  the  experimental  data.  The 
resultant  heat  transfer  rate,  Q,  is  then  substituted  into  equation  (5.9)  to  find  the 
overall  heat-transfer  coefficient,  U^. 


Q 

A^iLMTD) 


(5.12) 


where: 

Q  =  heat  transfer  rate  from  eq.  (5.11)  (W) 

Ao  =  effective  outside  surface  area  (m^) 

LMTD  =  log  mean  temperature  difference;  eq.  (5.10)  (K) 
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Since  U^,  A^,  and  A,  are  known  quantities,  this  leaves  only  two  unknowns 
in  equation  (5.8),  and  h„  the  outside  and  inside  heat  transfer  coefficients. 

Often  the  coolant  side  thermal  resistance  is  dominant,  and  inserts  such  as  those 
mentioned  in  section  IV  C  are  used  to  lower  the  inside  resistance.  This  allows  a 
more  accurate  computation  of  the  outside  heat-transfer  coefficient,  h,,,  when  using 
the  modified  Wilson  plot  technique  mentioned  in  section  V  B.  Vapor  side  heat 
transfer  may  also  be  enhanced  through  the  use  of  fins,  drainage  strips,  or  dropwise 
condensation  promoters. 

B.  MODIFIED  WILSON  PLOT  TECHNIQUE 

The  ideal  way  to  solve  for  h,,  and  h,  in  equation  (5.8)  is  through  the  use  of 
instrumented,  which  accurately  determine  a  mean  tube  wall  temperature.  The  inside 
and  outside  mean  tube  wall  temperatures  may  then  be  obtained  directly  by  assuming 
a  linear  temperature  profile  across  the  wall.  Since  the  vapor  temperature  and  mean 
coolant  temperature  are  known,  the  inside  and  outside  heat-transfer  coefficients  may 
then  be  calculated  directly  using  equation  (5.13). 

q=h^T  (513) 

where: 

q  =  heat  flux  (W/m^) 

h  =  heat-transfer  coefficient  (h;  for  inside,  h^  for  outside)  (W/m^  K) 

AT  =  temperature  difference  across  resistive  medium  (AT  =  T„j-T,^,u,puuide 
for  the  vapor  side,  and  AT  =  T„u.i„,ide  -  for  the  coolant  side) 

For  data  collection  on  a  large  number  of  tubes,  the  use  of  instrumented  tubes 
is  impractical  due  to  the  high  cost  and  difficulty  involved  in  manufacturing  so.  many 
tubes.  Therefore,  the  modified  Wilson  plot  technique  was  developed,  which  solves 
for  the  inside  and  outside  heat-transfer  coefficients  simultaneously  without  using  wall 
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thermocouples.  To  obtain  the  most  accurate  results  with  this  method,  it  is  necessary 
that  the  inside  and  outside  coefficients  be  relatively  equal  in  magnitude. 

The  modified  Wilson  plot  technique  requires  that  the  "form"  of  the  equation 
for  both  the  inside  and  outside  heat-transfer  coefficients  be  known.  The  Nusselt 
theory  and  Sieder-Tate  correlation  are  used  to  represent  the  "form"  of  the  outside 
and  inside  heat-transfer  coefficients  respectively.  The  Nusselt  theory  when  based  on 
q  can  be  represented  by: 


h 


O 


(5.14) 


where: 

ho  =  outside  heat-transfer  coefficient  (W/m^  K);  based  on  q 
a  =  dimensionless  coefficient 

=  thermal  conductivity  of  condensate  film  (W/m  K) 

Pf  =  condensate  film  density  (k^m^) 

g  =  gravitational  constant  (9.81  m/s^) 

hfg  =  specific  enthalpy  of  vaporization  (J/kg) 

=  dynamic  viscosity  of  condensate  film  (N  s/m^) 

Do  =  outside  tube  diameter  (m) 

q  =  heat  flux  based  on  outside  area  (Q/Ao)  (W/m^) 

The  Sieder-Tate  correlation  may  be  represented  by: 


h,=C,-^Re°*Pr''^ 

•D, 


\0.14 


=qQ 


(5.15) 


where: 

h,  =  inside  heat  transfer  coefficient  (W/m^  K) 

C,  =  Sieder-Tate  leading  coefficient 

ko  =  thermal  conductivity  of  cooling  water  (W/m-K) 

D,  =  inside  tube  diameter  (m) 

Re  =  Reynolds  number 

Pr  =  Prandtl  number 

/lo  =  dynamic  viscosity  of  cooling  water  at  bulk  temperature  (N  s/m^) 
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=  dynamic  viscosity  of  cooling  water  at  mean  inner  wall  temperature 
(N*s/m^) 

Substituting  equations  (5.14)  and  (S.IS)  into  equation  (5.8)  gives: 


(5.16) 


By  letting: 


-R^A 


Z 


(5.17) 


(5.18) 


m  =  J-  (5.19) 

f>  =  l  (5.20) 

a 

then  equation  (5.16)  reduces  to: 

Y=mX^b  (5-21) 

The  parameters  H  and  Z  are  both  temperature  dependent,  and  must  be 
determined  iteratively.  A  least-squares  fit  of  equation  (5.21)  is  used  to  determine 
C,  anda.  Once  C,  is  known,  h,  can  be  calculated  using  equation  (5.15).  With  h,  and 
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Uo  known,  the  value  of  h„  can  then  be  easily  determined  by  rearranging  equation 
(5.7)  as  in  equation  (5.22),  or  by  using  equation  (5.14)  with  the  value  of  a  known. 


(5.22) 


C.  INSTRUMENTED  TUBE  IMPROVEMENTS  FOR  DATA  REDUCTION 

In  previous  work  at  NPS,  the  standard  form  of  the  Sieder-Tate  equation  was 
used  with  a  Reynolds  number  exponent  of  0.8,  as  in  equation  (2.10).  One  of  the 
aims  of  this  thesis  was  to  use  an  instrumented  tube  to  directly  determine  the  inside 
and  outside  coefficients,  h,  and  h„  ,  and  then  use  the  data  to  determine  a  more 
"exact"  form  of  the  Sieder-Tate-type  equation  to  be  used  for  each  insert.  The 
coolant  side  correlations  mentioned  in  section  II  C  were  based  on  a  long,  straight 
entrance  length.  The  sharp  90°  bend  just  prior  to  the  test  section  tube  entrance 
undoubtedly  creates  entrance  effects  which  lead  to  discrepancies  between  our 
experimental  data  and  heat  transfer  behavior  predicted  by  the  Sieder-Tate 
correlation. 

Assuming  the  final  form  of  the  inside  heat-transfer  correlation  to  be  that  of 
equation  (2.10)  gives  the  following: 


(5.23) 


where: 

m  =  Reynolds  number  exponent  to  be  determined 
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Rearranging  equation  (5.23)  gives: 


Nu 


Pr 


(5.24) 


Taking  the  natural  log  of  equation  (5.24)  yields: 


In 


Nu 

f  IL  > 

0.14 

p^l/3 

Pc 

> 

kk  Ci+  m  hi  Re 


(5.25) 


Equation  (5.25)  is  in  the  form  of  a  linear  equation,  and  by  plotting  In 
(Nu/(Pr‘^(/i</m,)® ‘^))  versus  In  Re,  the  slope  and  intercept,  namely  the  Reynolds 
number  exponent  and  Sieder-Tate  coeifficient  may  be  determined  from  the 
instrumented  data. 

With  the  unknown  parameters  of  equation  (5.23)  determined,  the  new  inside 
heat-transfer  correlations  (one  for  each  insert)  could  be  used  in  the  data  reduction 
program  to  give  the  value  of  the  inside  heat-transfer  coefficient,  h„  directly  and  to 
provide  a  more  accurate  calculation  of  the  outside  heat-transfer  coefficient  h^. 

D.  ENHANCEMENT  RATIO 

Following  the  development  of  Van  Petten  [Ref.2],  and  Nusselt  theory, 
experimental  data  can  be  curve  fitted,  using  a  least-squares  analysis,  to  an  equation 
of  the  following  form: 
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(5.26) 


q  =  aAT" 


where: 


It 


=  heat  flux  based  on  outside  area  (Q/A„)  (W/m^ 
=  temperature  drop  across  the  condensate  nlm  (K) 

Substituting  equation  (5.12)  into  his  expression  yields: 


h  =  a  AT"'* 

O 


(5.27) 


From  Nusselt  theory  n  =  0.75,  therefore  the  enhancement  ratio,  based  on 
constant  temperature  drop  across  the  condensate  fllm,  can  be  expressed  as: 


€ 


Ar  " 


(5.28) 


where: 

=  enhancement  ratio  based  on  constant  temperature  drop  across  the 
condensate  film 

/  =  subscript  denoting  finned  tube 

s  =  subscript  denoting  smooth  tube 

hj,  =  outside  heat-transfer  coefficient  (W/m^  K) 

a  =  constant  of  proportionality  introauced  in  equation  (5.26) 


Also  for  constant  AT;  using  equation  (5.14): 


€ 


AT 


^  o, 

a/iT 

(5.29) 


If  the  heat  flux  is  kept  constant  the  values  of  Zj  and  Z,  remain  equal,  which 
results  in  equation  (5.30). 
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e  -  V|  - 
“  «  a.Z.  o. 

At  S  S  S 


(5.30) 


Combining  equations  (5.29  and  (5.30)  gives  the  relationsh^  between  and 
€  in  equation  (5.33). 


II 

s 

(5.31) 

«4^/j 

II 

Jvft 

(5.32) 

[“.J 

(5.33) 

Note  that  and  e  are  independent  of  q  and  AT. 
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VI.  RESULTS  AND  DISCUSSION 


A.  DROPWISE  CONDENSATION 

As  mentioned  in  section  IV  B,  a  dropwise  condensation  condition  was  obtained 
initially.  Data  were  taken  during  dropwise  conditions  at  vacuum  and  atmospheric 
pressure  with  the  instrumented  tube,  and  this  data  is  compared  to  filmwise  data  in 
Figure  4.  The  figure  shows  a  marked  contrast  between  filmwise  and  dropwise 
condensation  data.  The  dropwise  heat-transfer  enhancement,  compared  to  the 
filmwise  data,  varied  from  ~2  to  7  for  vacuum  conditions  and  ~9  to  10  for 
atmospheric  conditions.  Marto  et  al  [Ref.  23]  studied  the  use  of  organic  coatings  for 
the  promotion  of  dropwise  condensation  of  steam,  and  obtained  an  outside  heat- 
transfer  coefficient  of  ~55  kW/m^  K  for  a  Fluoroacrylic  coating  compared  to  ~30  to 
85  kW/m^  K  found  under  vacuum  conditions  in  this  study  (both  taken  at  P~  11  kPa 
and  vapor  velocity  ~2  m/s).  The  difficulty  in  accurately  measuring  the  temperature 
drop  across  the  condensate  film  for  the  dropwise  condition  is  illustrated  by  the  large 
amount  of  scatter  in  the  dropwise  data;  this  problem  is  caused  by  the  lack  of  a  stable 
film  and  the  intermittent  presence  of  drops  near  the  instrumented  tube  wall 
thermocouples.  The  exact  cause  of  the  dropwise  conditions  was  never  determined. 
Clearly,  some  organic  contamination  either  from  the  boiler  feed  water  or  fi’om  the 
gasket  material  was  depositing  on  the  test  tube.  Since  this  thesis  was  devoted  to 
filmwise  condensation,  great  efforts  were  made  to  clean  the  test  tube  and  prepare 
it  chemically  so  that  the  condensate  would  wet  the  surface. 
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Figure  4.  Comparison  of  Dropwise  and  Fifanwise  Condensation  Data  (Smooth 
Instrumented  Tube,  Heatex  Insert) 


48 


B.  INSTRUMENTED  TUBE  RESULTS 


After  resorting  to  the  oxide  coating  procedure  mentioned  previously  to  mitigate 
the  dropwise  contamination  problem,  a  series  of  runs  were  made  with  the 
instrumented  tube  fabricated  by  Poole  [Ref.  6]  with  the  Heatex  insert  installed. 
Using  a  mean  waU  temperature,  the  inside  and  outside  heat-transfer  coefficients 
could  be  evaluated  directly. 

Figure  5  shows  filmwise  condensation  data  taken  at  various  pressures  and  vapor 
velocities.  The  four  data  sets  for  increasing  vapor  velocity  (from  1  to  3.5  m/s)  and 
decreasing  vapor  pressure  (from  101  to  28  kPa)  show  the  effect  of  vapor  shear 
thinning  the  condensate  film,  giving  an  increase  in  the  outside  heat-transfer 
coefficient.  This  series  of  data  runs  was  taken  by  maintaining  the  heater  voltage  at 
175  volts  and  adjusting  coolant  flow  through  the  auxiliary  condenser  to  control 
pressure. 

The  four  data  sets  in  Figure  5  with  constant  vapor  velocity  (—1  m/s)  and 
decreasing  vapor  pressure  (from  101  to  28  kPa)  show  reduction  of  the  outside  heat- 
transfer  coefficient  with  decreasing  saturation  temperature  (due  to  decreasing 
saturation  pressure).  This  effect  is  thought  to  be  due  to  an  increase  in  condensate 
viscosity  at  lower  temperatures  which  tends  to  prevent  the  condensate  from  flowing 
around  and  draining  from  the  tube  as  easily  as  at  higher  temperatures.  The  resulting 
condensate  film  thickening  provides  an  additional  resistance  to  heat  transfer  thereby 
lowering  the  outside  heat-transfer  coefficient.  This  series  of  data  runs  was  taken  by 
both  adjusting  heater  power  and  auxiliary  condenser  coolant  flow  to  obtain  the  same 
pressure  conditions  as  above  with  a  constant  vapor  velocity  of  ~  1  m/s. 
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( 3.5  m/s ,  28  kPa ) 


(  >1  •  2vUi  /  M  )  OM 


Figure  5.  Effect  of  Pressure  and  Vapor  Velocity  on  the  Steam 

Heat-Transfer  Coefficient  (Smooth  Instrumented  Tube, 
Heatex  Insert) 
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Figure  6  shows  how  increased  pressure,  with  vapor  velocity  held  constant,  gives 
a  larger  increase  in  the  outside  heat-transfer  coefficient  over  the  predicted  Nusselt 
value.  This  discrepancy  is  also  thought  to  be  caused  by  the  viscosity  effect  since  the 
Nusselt  treatment  seemingly  takes  into  account  the  other  possible  causes.  Also 
depicted  more  clearly  is  the  dramatic  enhancement  due  to  vapor  shear  effects  for 
increased  vapor  velocity  at  constant  pressure. 

As  mentioned  previously,  the  use  of  inserts  allows  more  effective  mixing  of 
coolant  and  facilitates  greater  accuracy  in  the  calculation  of  the  outside  heat-transfer 
coefficient.  Figure  7  shows  the  outside  heat-transfer  coefficient  determined  from  the 
use  of  the  two  inserts  together  with  no  insert.  The  results  for  the  wire  wrap  and 
Heatex  inserts  are  closely  grouped,  whereas  the  no  insert  data  shows  much  greater 
scatter.  This  seems  to  indicate  that  the  use  of  an  effective  insert  does  indeed 
enhance  accuracy  for  the  instrumented  tube  data. 

The  mean  temperature  difference  across  the  condensate  film  for  the  no  insert 
case  has  a  significantly  lower  value  than  either  of  the  insert  cases;  with  the  outside 
and  tube  wall  conditions  unchanged,  this  indicates  a  higher  inside  resistance  for  the 
no  insert  case.  To  illustrate  this  point  a  mid-range  instrumented  tube  data  point  at 
the  same  conditions  for  each  insert  case  is  shown  in  Table  3  (P„,~101  kPa, 
Vvapor~  11  ni/s,  V5j„^„,~2.75  m/s).  It  shows  that  both  inserts  have  a  comparable  effect 
(slightly  better  with  Heatex)  and  roughly  provide  a  factor  of  two  enhancement  over 
the  no  insert  case  in  the  inside  heat  transfer  coefficient.  Due  to  the  increased  inside 
resistance  for  the  no  insert  case,  the  heat  flux  shows  a  decrease  of  about  15%. 
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Figure  6.  Con^arison  of  Experimental  Results  with  Nusselt  Theory 
for  Varying  Pressures  and  Vapor  Velocities  (Smooth 
Instrumented  tube,  Heatex  Insert) 
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Figure  7.  Comparison  of  Steam  Heat-Transfer  Coefficients  at  Atmo^heric 

Conditions  (P„,~101  kPa,  m/s)  for  Three  Insert  Conditions 

(Smooth  Instrumented  Tube) 
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Table  3.  PERFORMANCE  COMPARISON  OF  TUBE  INSERTS 


AT  ATMOSPHERIC  CONDITIONS 


Heatex  Insert 

Wrap  Inser,^ 

No  Insert 

h,  (kW/m^  K) 

30.85 

27.63 

15.31 

q  (kW/m^) 

481.4 

477.7 

409.6 

h„  (kW/m^  K) 

9.86 

9.82 

10.75 

AT(K) 

48.82 

48.65 

38.10 

Figure  8  shows  the  wall  temperature  profiles  for  the  data  points  listed  in  Table 
3  along  with  a  vacuum  data  set  (P„,~28  kPa,  m/s,  Vj^,,„,~2.75  m/s)  for 

comparison,  where  0  degrees  is  at  the  top  dead  center  position  of  the  tube.  The 
shape  of  the  temperature  profiles  shows  the  effect  of  condensate  fihn  thickening 
toward  the  bottom  of  the  tube.  The  higher  resistance  through  a  thicker  condensate 
film  results  in  a  lower  tube  wall  temperature  toward  the  bottom  of  the  tube  as 
shown. 

As  expected,  the  wire  wrap  insert  and  Heatex  insert  temperature  profiles  at 
atmospheric  pressure  are  very  similar.  The  no  insert  profile  has  a  shape  similar  to 
the  wire  wrap  and  Heatex  profiles,  but  the  mean  wall  temperature  is  about  11  K 
higher  and  shows  the  effect  of  the  increased  inside  resistance  to  heat  transfer.  The 
vacuum  run  temperature  profile  shows  the  effect  of  a  lower  temperature  gradient 
between  the  steam  and  coolant;  the  lower  heat-transfer  potential  results  in  lower 
heat  fluxes  and  a  flatter  temperature  profile. 

Figure  9  shows  the  comparison  of  the  instrumented  data  with  the  predictions 
of  Nusselt,  Fujii,  and  Shekriladze  and  Gomelauri  covered  in  section  II  B.  The  data 
depicted  ranges  from  P„,  =  28  kPa  and  =  3.5  m/s  to  P„,  =  101  kPa  and 
=  1.1  m/s.  The  data  seems  to  follow  the  Shekriladze  and  Gomelauri  prediction  the 
closest,  but  is  also  very  near  the  Fujii  prediction. 
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Figure  8.  Horizontal  Tube  Wall  Temperature  Profiles  (Smooth 
Instrumented  Tube,  m/s) 
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Figure  9.  Comparison  of  Experimental  Results  with  the  Predictions 
of  Nusselt,  Fujii,  and  Shekriladze-Gomelauri. 
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C.  INSIDE  HEAT-TRANSFER  CORRELATIONS  FROM  INSTRUMENTED 


TUBE  RESULTS 

The  instrumented  tube  provides  a  method  of  determining  an  inside  heat- 
transfer  coefficient  for  each  insert  from  direct  wall  temperature  measurements. 
Wanniarachchi  et  al  [Ref.  8],  in  a  similar  effort  to  resolve  the  differences  between 
the  Sieder-Tate  correlation  and  his  experimental  results,  developed  a  correlation 
based  on  the  Sieder-Tate  expression.  Using  an  intercept  form,  with  data  taken  on 
the  same  apparatus  as  used  in  this  study,  and  a  least-squares  fitted  leading 
coefficient,  the  correlation  took  the  following  form: 


=0.064 


+26.4 


(6.1) 


However,  Rouk  [Ref.  4],  using  an  optimization  technique,  showed  that  the  value  of 
the  intercept  had  little  effect  on  the  results,  and  it  was  the  accurate  determination 
of  the  Reynolds  exponent  that  was  more  critical. 

With  a  view  to  finding  the  appropriate  Reynolds  exponent  and  leading 
coefficient  for  equation  (5.23),  the  data  for  no  insert,  wire  wrap,  and  Heatex  insert 
were  plotted  as  In  Re  versus  In  (Nu/Pr‘^(/Ltyju*)®  ”)  as  explained  in  section  V  C.  The 
plotted  data  are  shown  in  Figures  10,  11  and  12  for  no  insert,  wire  wrap  insert,  and 
Heatex  insert  respectively.  A  line  of  best  fit  (typically  with  a  regression  coefficient 
of  0.99)  was  used  to  obtain  the  value  of  the  intercept.  Figure  13  shows  the  plot  of 
aU  three  cases  on  the  same  graph  for  comparison.  Note  that  the  insert  data  is 
closely  grouped  on  the  upper  regions  of  the  graph  when  compared  to  the  no  insert 
data.  The  increase  in  the  Nusselt  number  again  indicates  more  efficient  inside  heat 
transfer  for  the  insert  vice  the  no  insert  case.  The  Reynolds  exponent,  or  slope. 
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Figure  10.  Log-Log  Plot  of  Re  versus  Nu/Pr^^^/i^)®  ”  for  No 
Insert  (Smooth  Instrumented  Tube) 
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Figure  11.  Log-Ijos  Plot  of  Re  versus  Nu/Pr*^(ay]^.*)“  ”  for  Wire  Wrap 
Insert  Smooth  Instrumented  Tube) 
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Figure  12.  Log-Log  Plot  of  Re  versus  Nu/Pr‘^(^a/im,)® for  Heatex 
Insert  (Smooth  Instrumented  Tube) 
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differs  from  the  Sieder-Tate-type  equation  of  0.8  in  all  three  cases.  The  following 
derived  correlations  apply  specifically  onty  to  the  medium  tube,  but  should  be 
applicable  to  any  tube  with  the  same  inside  diameter. 

The  no  insert  inside  heat-transfer  correlation  had  the  form: 


Ml  =0.013  R£^*^Pt 


(6.2) 


Equation  (6.2)  was  used  to  reprocess  some  current  smooth  tube  runs,  with  no  insert, 
as  well  as  those  of  Guttendorf  [Ref.  3]  and  Van  Petten  [Ref.  2]  to  check  their  values 
for  the  smooth  tube  Nusselt  coefficient  a. 

The  wire  wrap  insert  inside  heat-transfer  correlation  had  the  form: 


Aftt  =0.052 


(6.3) 


Equation  (6.3)  was  used  to  reprocess  previous  data  taken  on  the  same  apparatus  by 
Guttendorf  and  Van  Petten,  who  used  the  wire  wrap  insert  for  their  finned  and 
smooth  tube  experiments,  with  a  view  to  checking  their  results  with  this  new 
correlation.  Table  4  shows  a  comparison  between  the  new  wire  wrap  insert 
correlation  (eq.  6.3),  and  Wanniarachchi’s  correlation  (eq.  6.1)  both  with  and  without 
the  intercept  value  included.  This  comparison  shows  that  equation  6.1  is  actually 
more  accurate  without  the  intercept  value  included  when  compared  to  the  results  of 
the  new  wire  wrap  insert  correlation. 
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Table  4.  COMPARISON  OF  EQUATIONS  (6.1)  AND  (6.3) 

POP  PrW/,.  4 


FOR  Pr^(M7A.)®  ‘"=1.4  (held  cons 

tant  for  companion) 

Reynolds  Number 

20,000 

30,000 

40,000 

50,000 

Nusselt  Number 

Nu =0.064  Re“®Pr‘^(Aty/i.)‘“^-»-26.4 

273.7 

368.4 

456.9 

541.0 

Nusselt  Number 

Nu= 0.064  Re^^Pr'^CLtyA,^)®  ** 

247.3 

342.0 

430.5 

514.6 

Nusselt  Number 

Nu=0.052  Re^^Pr^^fA^M.)®” 

244.9 

341.5 

432.3 

519.1 

The  Heatex  insert  inside  heat- transfer  correlation  had  the  form; 


Nu=0.12  Re^^Pr^ 


\0.14 


V^w/ 


(6.4) 


Equation  (6.4)  was  used  to  reprocess  all  Heatex  data. 

Memory  [Ref.  12]  conducted  condensation  experiments  on  a  different  apparatus 
with  a  smooth  instrumented  tube.  He  also  determined  the  Reynolds  number 
exponent  for  no  insert,  wire  wrap  insert  (made  locally  and  somewhat  different  from 
the  one  used  in  this  study),  and  Heatex  insert.  The  exponents  he  obtained  are 
reported  in  Table  5.  The  Heatex  insert  exponent  of  0.68  compares  very  favorably 
with  the  value  of  0.69  found  in  this  study.  The  wire  wrap  exponent  of  0.73  was  well 
below  the  value  of  0.82,  most  likely  due  to  differences  in  the  insert.  The  no  insert 
case  gave  a  value  of  0.85  and  compared  well  to  the  value  of  0.89  found  in  this  study. 
The  difference  in  the  no  insert  exponent  is  most  likely  due  to  the  difference  in  the 
tube  entrance  region;  Memory  used  a  long  run  of  straight  pipe  for  the  tube  entrance, 
whereas  the  present  work  had  a  sharp  bend  just  prior  to  the  condenser  tube. 

Rouk  [Ref.  4]  used  the  instrumented  smooth  tube  data  of  Georgiadis  [Ref.  5] 
to  the  find  the  appropriate  value  of  the  Reynolds  number  exponent  for  the  no  insert 
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and  wire  wrap  cases.  His  results  compare  quite  well  with  this  study  and  are  also 
given  in  Table  5. 

ANL  (Argonne  National  Laboratoiy)  [Ref.  21]  conducted  an  assessment  of 
heat-transfer  correlations  for  turbulent  pipe  flow  with  water  to  determine  the  best 
correlation(s)  on  which  to  base  their  design  of  Ocean  Thermal  Energy  Conversion 
(OTEC)  heat  exchangers.  ANL  used  two  shell-and-tube  heat  exchangers,  with  no 
inserts,  for  analysis  and  reported  the  following: 

1.  The  Dittus-Boelter  (eq  2.7)  and  Sieder-Tate  (eq  2.9)  correlations  under¬ 
predicted  the  data  by  5%  to  15%  and  were  considered  too  conservative  for 
design. 

2.  Overall,  the  "best"  correlations  were  found  to  be  Petukhov-Popov  (eq  6.5)  and 
Sleicher-Rouse  (eq  6.6),  both  of  which  showed  excellent  agreement  (±  5%) 
with  the  experimental  data  (at  Pr=6.0  and  Pr=11.6). 


where; 


€ 

K, 


Nu  = 


(e/S)RePr 


K^+K^(elSy'^(Pr^-l) 

(valid  for  0.5  <  Pr  <  2000  and  10'  <  Re  <  5x10*) 


^1.82^1og,o  Re  -  1.64)-^ 
11.7  +  1.8  Pr  '^ 


Nw  =  5 +0.015  Re/Pri 

(valid  for  0.1  <  Pr  <  1(P  and  10'  <  Re  <  10*) 

where; 

a  =  0.88  -  0.24/a-l-PrJ 
b  =1/3-1-  0.5e-**‘^ 


(6.5) 


(6.6) 
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3.  The  most  accurate  correlations  (i.e.  Petukhov-Popov  and  Sleicher-Rouse) 
seem  to  yield  effective  Reynolds  exponents  in  the  neighborhood  of  0.85 
(uncertainty  range:  m  =  0.82  to  0.88). 

4.  The  potential  sources  of  uncertainty  in  the  Wilson  procedure  included 
waterside  flow  maldistribution,  entrance  effects,  experimental  error  in  U  ,  and 
the  uncertainty  in  the  Reynolds  number  exponent.  Of  these,  they  concluded 
that  the  uncertainty  in  the  Reynolds  number  exponent  was,  by  far,  the  most 
significant.  In  fact,  the  results  of  the  Wilson  procedure  were  found  to  be 
highly  sensitive  to  the  value  of  the  Reynolds  number  exponent. 


Table  5.  COMPARISON  OF  REYNOLDS  NUMBER  EXPONENTS  FOR 


SIEDER-TATE-TYPE  CORRELATIONS 


Experimental 

Data 

Rouk 

Memory 

ANL 

No  Insert 

0.89 

0.90 

0.85 

0.85 

Wire  Wrap 

0.82 

0.78 

0.73 

— 

Heatex 

0.69 

— 

0.68 

— 

D.  ANALYSIS  OF  SMOOTH  TUBE  RESULTS 

When  using  the  modified  Wilson  plot  technique  to  reprocess  data  files,  the 
solution  option  can  be  specified  to  use  either  the  stored  value  of  the  Sieder-Tate 
coefficient  (for  direct  computation  of  h^)  or  let  the  coefficient  value  "float",  which 
allows  the  program  to  calculate  its  own  value  of  the  coefficient.  In  order  to 
determine  which  method  was  most  accurate,  the  instrumented  data  files  were 
reprocessed  using  each  method  and  then  compared  with  the  values  of  the  heat- 
transfer  coefficient  which  were  obtained  by  direct  measurement  of  the  tube  wall 
temperature.  A  high,  medium,  and  low  coolant  flow  rate  was  chosen  from  each  run 
to  facilitate  the  comparison.  The  results  were  tabulated  and  are  shown  in  Appendix 
D;  it  can  be  seen  that  the  fixed  coefficient  method  yielded  the  more  accurate  results 
at  least  75%  of  the  time.  The  mean  error  of  the  fixed  method  was  ±2.0%,  and  that 
of  the  floating  method  was  ±5.4%.  The  error  for  the  lowest  coolant  flow  rate  (Re 
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<20,000)  was  noticeably  higher  than  for  higher  coolant  flow  rates  for  both  methods. 
The  choice  of  using  the  fixed  coefficient  method  represents  a  departure  from  the 
practice  of  previous  researchers  on  this  apparatus  who  exclusively  used  the  floating 
coefficient  method. 

Prior  to  the  instrumented  tube  runs  reported  in  Figure  5  (from  which  the  new 
correlations  were  empiricaUy  derived)  a  series  of  data  runs  were  made  using  a  plain 
smooth  uninstrumented  tube  (S02)  of  the  same  dimensions  using  the  Heatex  insert. 
The  plain  smooth  tube  data  was  then  reprocessed  using  the  new  correlation  for  the 
Heatex  insert  and  the  results  are  plotted  in  Figure  14.  These  data  sets  were  taken 
at  the  same  conditions  as  the  instrumented  data  of  Figure  5  except  for  the  set  at  the 
highest  vapor  velocity  of  6.2  m/s  vice  3.5  m/s  for  the  instrumented  tube.  Similar 
effects  of  vapor  shear  and  vapor  pressure,  as  mentioned  previously  for  Figure  5,  are 
clearly  seen,  and  again  illustrate  the  vapor  shear  effect  on  the  outside  heat>transfer 
coefficient. 

With  the  exception  of  the  two  data  runs  at  high  vapor  velocity,  the  data  fi’om 
Figures  5  and  14  are  shown  together  in  Figure  15.  The  close  agreement  of  the 
reprocessed  plain  smooth  tube  data  with  the  instrumented  smooth  tube  data  aUows 
a  high  degree  of  confidence  in  the  accuracy  of  the  new  correlations  and  the  choice 
of  the  fixed  coefficient  method. 

To  provide  a  baseline  from  which  to  evaluate  finned  tube  performance  it  was 
necessary  to  obtain  the  smooth  tube  Nusselt  coefficient,  a,  for  the  specific  conditions 
under  which  the  comparison  was  to  be  made.  The  condition  chosen  was  atmospheric 
pressure  (101  kPa)  and  a  vapor  velocity  of  ~1  m/s. 

For  8  complete  sets  of  data  the  average  value  of  a,  using  the  fixed  method,  was 
found  to  be  0.876.  The  average  value  for  each  data  set  was  found  by  taking  the 
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Figure  14.  Effect  of  Pressure  and  Vapor  Velocity  on  the  Steam 
Heat-Transfer  Coefficient  (Non-instrumented  Smooth 
Tube,  Heatex  Insert) 
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Figure  15.  Comparison  of  Instrumented  Smooth  Tube  Results  with  Non- 
Instrumented  Smooth  Tube  Data  After  Reprocessmg  with  the 
New  Heatex  Insert  Inside  Heat-Transfer  Correlation. 
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measured  value  of  for  each  data  point,  dividing  it  by  the  Nusselt  theory  prediction 
of  ho,  and  then  multiplying  by  0.728  (the  Nusselt  coefficient);  the  average  of  all  the 
data  points  in  the  set  was  then  taken.  Interestingly  the  Wilson  plot  floating 
coefficient  method  gave  a  value  of  0.835,  somewhat  lower  than  the  average  value. 
Originally  it  was  thought  that  this  discrepancy  might  be  due  to  "outlier"  data  points 
(high  or  low  coolant  velocity)  in  each  set.  However,  removing  the  highest  or  lowest 
coolant  flow  rates  within  a  set  had  little  or  no  effect  on  the  Wilson  plot  result.  The 
reason  for  the  discrepancy  is  still  not  known  and  merits  future  study. 

The  value  of  a  was  calculated  for  several  other  flow  conditions;  these  are 
shown  in  Table  6.  The  trend  of  the  readings,  like  that  of  Figure  6,  shows  that  vapor 
velocity  has  a  much  greater  effect  than  pressure  on  the  value  of  a,  as  expected. 


Table  6.  SMOOTH  TUBE  a  SUMMARY;  EFFECT  OF  PRESSURE  AND 
VAPOR  VELOCITY 


File  Name 

P(kPa) 

Vapor  Velocity 

(m/s) 

a 

FlMAVSHl 

28 

3.5 

1.015 

FIMAVSH2 

41 

2.5 

.985 

FIMAVSH3 

68 

1.5 

.930 

FIMAVSH4 

101 

1.1 

.866 

FIMAVSH7 

69 

1.0 

.836 

FIMAVSH6 

41 

1.0 

.818 

FIMAVSH5 

28 

1.0 

.786 

E.  ANALYSIS  OF  nNNED  TUBE  RESULTS 

With  an  accurate  value  of  a,  and  the  newly  determined  inside  heat-transfer 
correlations,  the  medium  family  finned  tube  data  of  Van  Petten  [Ref.  2]  was 
evaluated.  Figure  16  shows  the  data  Van  Petten  reported  in  his  thesis;  it  also  shows 
his  data  after  being  reprocessed  using  tne  new  wire  wrap  insert  correlation  with 
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Figure  16.  Comparisou  of  the  Steam  Heat-Transfer  Enhancement  Data  of  Van 
Petten,  for  the  Medium  Finned  Tube  Family,  Using  the  Modified 
Wilson  Plot  and  New  Wire  Wrap  Insert  Inside  Correlation. 
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fin  spacing  (  mm 


both  the  fixed  coefficient  and  Wilson  plot  floating  coefficient  methods.  Since  Van 
Petten  used  the  Wilson  plot  method,  it  is  not  surprising  that  the  original  thesis  data 
and  new  Wilson  plot  floating  coefficient  data  are  comparable  since  the  Reynolds 
exponent  only  varied  from  0.8  to  0.82.  The  fixed  coefficient  method  enhancement 
is  substantially  higher  than  the  Wilson  plot  results. 

Since  the  assertion  is  that  the  fixed  coefficient  method  is  more  accurate  than 
the  Wilson  plot  method,  then  the  conclusion  must  be  that  the  enhancement  for  this 
set  of  finned  tubes  is  actually  higher  than  previously  reported. 

During  this  study,  limited  medium  finned  tube  experiments  were  conducted  for 
purposes  of  comparison.  Figure  17  shows  the  comparison  between  this  data  and  the 
newly  reprocessed  data  of  Van  Petten  (using  the  fixed  coefficient  method)  and  shows 
reasonable  agreement.  To  more  clearly  illustrate  this  point,  the  data  taken  on  the 
2.0  mm  fin  spacing  tube  has  been  given  in  more  detail  in  Figure  18.  Excellent 
agreement  is  seen  between  the  experimental  results  of  this  study  and  that  of  Van 
Petten  using  the  known  inside  heat-transfer  correlation  with  the  fixed  coefficient 
method.  Again,  as  shown  in  Figure  16,  the  Wilson  plot  prediction  is  significantly 
below  the  fixed  coefficient  results. 
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Figure  17.  Comparison  of  the  Steam  Heat-Transfer  Enhancement  Data 

of  Van  Petten  and  Swensen  for  the  Medium  Finned  Tube  Famify. 


72 


spacing  (  mm 


o  ^ 

CO 


CB 

I- 


fO 

o  w 

CM 


O 

ir> 


o 

o 

o 

lO 


o 

o 

o 

o 

■M- 


o 

o 

o 

lO 

CO 


(  >i  *  3vUi  /  M  )  OM 


Figure  18.  Comparison  of  the  Steam  Heat  Transfer  Data  of  Van  Petten  and 
Swensen  for  the  2.0  mm  Fin  Spacing  Medium  Tube. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 


1.  The  inside  heat- transfer  correlation  is  highly  sensitive  to  the  Reynolds  number 
exponent. 

2.  Each  insert  condition  must  be  analyzed  separately  to  determine  the 
appropriate  "form"  of  the  inside  heat-transfer  correlation. 

3.  Calculations  based  on  a  known  inside  heat-transfer  correlation  are  more 
accurate  than  modified  Wilson  plot  results. 

4.  Armed  with  accurate  inside  heat-transfer  correlations,  previous  data  may  be 
reprocessed  to  give  more  accurate  results. 

5.  The  source  of  contamination  in  the  test  apparatus,  which  has  caused  a 
dropwise  condensation  problem  for  a  number  of  years,  is  most  probabfy  due 
to  a  contaminated  distilled  water  source. 


B.  RECOMMENDATIONS 


1.  Reprocess  aU  previous  medium  and  large  diameter  finned  tube  data  using  the 
fixed  coefficient  method  to  obtain  more  accurate  results. 

2.  Continue  with  construction  of  smooth  instrumented  tubes  of  different 
diameters  (i.e.  small,  medium,  and  large)  to  confirm  die  medium  tube  results 
and  develop  correlations  specifically  wt  the  small  and  large  diameter  tubes. 

3.  Construct  one  representative  instrumented  finned  tube  to  test  the  validity  of 
applying  instrumented  smooth  tube  results  to  finned  tube  data. 

4.  Test  representative  water  samples  that  have  been  collected  from  both  the 
distiller  and  boiler  to  confirm  the  presence  of  impurities  and  validate  their 
origin. 

5.  Replace  current  distiller  with  a  deionized  pure  water  source  (either 
commercial  purchase  of  water  or  new  distilling  apparatus). 
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APPENDIX  A.  PHYSICAL  AND  THERMODYNAMIC  PROPERTIES  OF 
WATER 

The  physical  and  thermodynamic  properties  of  water  were  based 
following  equations: 

Heat  Capacity  (J/kg-K): 

Cp  =  4.2 1 1  -  r*  [2.268;t  lO*^  -  T  *  (4.424 jc  lO^  +  2.714x  10’  *  7)] 
where:  T  =  temperature  (celsius) 

Dynamic  viscosity,  (kg/m  s): 

Thermal  conductivity  (W/m-K): 

k  =  -0.9225 +JC  *(2.8395 -x  *(1.8007  -x  *(0.5258  -0.0734  *x))) 
where:  x  =  (T+273.15)/273.15 

Density  (kg/m^): 

p  =999.5295  +  T*  (0.0127-7*  (5.4825xi0-’-7*  1.2341x70’)) 


on  the 


(A.1) 


(A.2) 


(A.3) 


(A.4) 
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Latent  heat  of  vaporization  (J/kg): 


=2477200  -2450  *(7- 10) 


Fluid  Enthalpy  (J/kg): 


T*  (4.2038  -  T*  (5.8813jc/0‘^-r*4.5516xi0'^) 


Saturation  pressure  (Pa): 


P^= 22120000  *77- 

Jw 


where: 


Pr 

=  e®' 

Br 

=  SUM  /rr,  (1+  4.1671  (1-T,)  + 

20.9751  (1-7,)^) 

((l.Ox  10'0( 

1-T,)^  +  6)] 

SUM 

=  (-7.6912)(1-T, 

)- 26.0802(1-7,)^- 

168.1707  (1-7,)^ 

118.9646  (1- 

ty 

T, 

=  (T+273.15)  / 

647.3 

(A.5) 

(A.6) 

(A.7) 

•  (1-T,)/ 

+  64.2329(1-7,)'- 
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APPENDIX  B.  SYSTEM  CALIBRATIONS  AND  CORRECTIONS 
B.1  Thermocouple  and  Quartz  Thermometer  Calibration 

Several  different  thermocouples  and  the  quartz  crystal  thermometer  (HP 
2804 A,  Ser.  No.  2244AD1192)  were  calibrated  against  a  platinum  resistance  probe 
in  a  mixed  isothermal  ethylene  glycol  bath  with  a  Rosemont  Galvanometer  model 
920A  commutating  bridge  (Ser.  no.  013494)  from  23  to  26  September,  1991. 

NOMENCLATURE 

T1  Quartz  Thermometer  Measuring  Probe  T1  (Ser.  No.  2120A-00707,  with  dial 

setting  481) 

T2  Quartz  Thermometer  Measuring  Probe  T2  (Ser.  No.  2120A-60459  with  dial 

setting  510) 

T1-T2  The  T1-T2  reading  on  the  quartz  thermometer 

T-55  Large  diameter  metal  sheath  thermocouple  (diameter  =  0.040") 

T-56  Small  diameter  metal  sheath  thermocouple  (diameter  =  0.020") 

T-57  Old  Thermocouple  (taken  from  rig  during  disassembly  set  to  HP  3497A 
internal  zero);  (Type:  Omega,  TT-T~30,  Lot#  HCP093HC0306) 

T-58  New  1  Thermocouple  (set  to  HP  3497A  internal  zero);  (Type:  Omega,  7T- 
T-30  SLE,  Lot#  OCP1453PTCC01473P) 

T-59  New  2  Thermocouple  (referenced  to  ice  bath  zero) 

T-60  10  Junction  Thermopile  #1  (referenced  to  ice  bath  zero) 

T-61  10  Junction  Thermopile  #2  (referenced  to  ice  bath  zero) 

T-62  10  Junction  Thermopile  #3  (referenced  to  ice  bath  zero) 

The  calibration  was  performed  by  taking  the  bath  temperature  up  from  290  K 
to  393  K  by  5  K  increments  then  back  down  to  check  for  hysteresis;  no  hysteresis  was 
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observed.  The  data  was  fitted  to  a  fifth-order  polynomial  (regression  coefficient  = 
1.000  for  each  polynomial  fit)  in  each  case  with  the  foUowing  results: 


755  =  273.15+(2.5943e-2)K-(7.2671-7)K2 

+  (3.2941-11)^3 -(9.7119-16)K*  (Bl) 

+  (9.7121 -20) 

756  =  273.15+(2.5878e-2)l'-(5.9853e-7)F2 

-  (3.1242e-ll)F3+(1.3275e-14)F‘  (B-2) 

-  (1.0188e-18)K3 

757  =  273.15+(2.5923e-2)F-(7.3933e-7)K3 

+  (2.8625e-ll)F3+(1.9717e-15)F^  (B-3) 

-  (2.2486e-19)F3 

T58  =  273.15+(2.5931e-2)F-(7.5323e-7)F^ 

+  (4.0567e-ll)F3-(1.2791e-15)F*  (B.4) 

+  (6.4402e-20)F’ 

759  =  273.15+(2.547U-2)F-(3.7621e-7)F2 

-  (1.0105«-10)F3+(2.3928e-14)F^  (B.5) 

-  (1.6440e-10)F3 

T60  =  273.17+(2.5571e-2)F-(1.9980«-7)K2 

-  (1.6385e-10)F3+(2.6164e-14)F^  (B.6) 

-  (1.0295e-18)F3 

T61  =  273.15+(2.6119e-2)F-(9.0449e-7)F^ 

+  (1.1214«-10)F3-(1.5623e-14)V^  (B.7) 

+  (1.0646^-18)1^ 
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T62  =  273.15+(2.5996e-2)K-(7.4405e-7)K2 

+  (2.4733c-ll)F^+(3.3236e-15)V*  (®-8) 

-  (3.7460tf-19)I^ 

For  T55,  T56,  T57,  T59\  V=  Voltage  in  microvolts  (eg.  0.010  volts  =  10,000 
microvolts). 

For  TOO,  T61,  T62-,  V=  Voltage  in  microvolts/10  (eg.  100,000  microvolts  10 

=  10,000). 

The  HP  2804A  Quartz  Crystal  Thermometer  was  also  cahbrated  on  26 
September,  1991  with  the  results  summarized  in  Table  B.l. 


Table  B.l.  QUARTZ  THERMOMETER  CALIBRATION  DATA 


Reference 

T1 

Error 
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Error 

Temperature 
(deg  C) 

(deg  C) 

(deg  C) 

(deg  C) 

(dec  C) 

17.979 

18.003 

0.024 

18.001 

0.022 

18.977 

18.997 

0.020 

18.996 

0.021 

19.984 

20.001 

0.017 

20.002 

0.018 

20.980 

21.007 

0.027 

21.008 

0.028 

22.936 

22.959 

0.023 

22.961 

0.025 

25.014 

25.032 

0.018 

25.037 

0.023 

30.140 

30.153 

0.013 

30.164 

0.024 

34.977 

34.983 

0.006 

34.998 

0.021 

The  mean  error  for  the  Quartz  Thermometer  as  shown  in  the  Table  B.l  is 
~+0.02  deg  C.  However,  it  is  well  noted  that  T1-T2,  the  critical  measurement, 
indicates  an  apparent  error  of  less  than  +0.005  deg  C  when  the  temperature  reading 
is  near  or  below  25  deg  C;  T1  and  T2  tracking  well  together  lowers  the  error 
estimate. 
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B.2  Flow  Meter  Calibration 

The  flow  meter  calibration  for  coolant  flow  through  the  single  horizontal  tube 
was  completed  on  Oct  28,  1991  using  a  stop  watch,  portable  tank,  and  a  Toledo 
model  31-0851  IV,  Se.  No.  1326  scale  with  1/10  pound  graduations.  The  following 
relation  was  obtained  via  linear  regression: 


m  =  6.7409F  +  13.027 


B.9 


where: 

m  =  mass  flow  rate  /grams  per  second) 

F  =  flow  meter  reading  (eg.  10%=^  10) 

The  applicable  range  of  the  calibration  was  10%  to  95%.  The  water 
temperature  on  the  day  of  the  calibration  was  17.5°C  (290.6  K)  and  water  density  was 
998.5  kg/m^.  The  data  is  shown  in  Figure  B.l 

B.3  Mass  Flow  Rate  Correction 

The  inlet  water  temperature  from  the  cooling  water  sump  varies  anywhere  from 
15°C  to  25°C  depending  on  environmental  conditions.  To  account  for  these 
temperature  variations  the  following  function  was  used  to  calculate  a  correction 
factor  for  viscosity  variation  with  temperature  [Ref.  5]. 

C/=1.0365  -(1.9644£-3)7In  +(5.2500e-6)7in2 

where: 

Cf  =  mass  flow  rate  correction  factor 
Tm  =  inlet  temperature  (celsius) 
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Figure  B.l  Horizontal  Tube  Coolant  Flowmeter  Calibration  Chart 
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The  value  of  Cf  for  the  flow  meter  calibration  (Tin=11.5°C)  was  1.0037. 
Therefore,  the  actual  mass  flow  rate  was  calculated  using  the  foUowing  equation; 


Cf 

1.0037 


(B.ll) 


where: 

mgc  =  corrected  mass  flow  rate  (grams  per  second) 
n^caic=  computed  mass  flow  rate  (eq.  B.9)  (grams  per  second) 


B.4  Pressure  Transducer  Calibration 

Three  methods  of  pressure  measurement  were  available  on  the  apparatus: 

1.  Direct  pressure  reading  off  the  Heise  solid  front  -  CM-104119  pressure  gauge, 
(range  0-15  psia). 

2.  Converted  voltage  readings  from  the  Setra,  model  204,  Ser.  no.  63982 
pressure  transducer  (range  0-14.7  psia;  0-5  volts;  5V~0  psia). 

3.  Steam  saturation  temperature  measurement  with  the  apparatus  producing 
steam  at  steady  state.  The  steam  saturation  temperature/pressure  relation 
was  utilized  via  standard  steam  tables. 

The  pressure  transducer  was  calibrated  versus  the  vapor  temperature  probe 
reading  on  12  December  1991.  Equation  B.12  gives  the  desired  relationship.  The 
data  is  shown  in  Figure  B.2. 


P  =  -2.9360K+ 14.7827 


(B.12) 


where: 

P  =  pressure  (psia) 

V  =  pressure  transducer  voltage  reading  (volts) 


B.5  FRICTION  TEMPERATURE  CORRECTION 

As  coolant  flows  through  the  tube  there  is  a  bulk  temperature  rise  due  to 
frictional  heating,  which  is  highly  dependent  on  fluid  velocity.  Although  small,  this 
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Figure  B.2 

Pressure  Transducer  Calibration  Chart 
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can  have  a  significant  effect  on  the  calculated  overall  heat-transfer  coefficient 
Measurements  were  made  for  no  insert,  Heatex,  wire  wrap  and  twisted  tape  inserts 
as  shown  in  Figure  B.3,  on  5  December  1991.  Each  data  set  was  curve  fitted  to  a 
third  order  polynomial  which  is  depicted  in  Table  B.2.  Each  respective  polynomial 
corrects  the  temperature  rise  measurement  for  the  heating  due  to  the  particular  type 
of  insert  used. 


Table  B.2  FRICTION  TEMPERATURE  RISE  POLYNOMIALS 


Insert  Tvne 

Polynomial 

None 

Trise  =  -1.960x10  ^V’+9.349xl0^V2-i-L749xl0^V-2.728xl0^ 

Wire  Wrap 

Trise  =  8.160xl0^V"+ 1.4512x10  V+2.745xl0*V-3.991xl0-^ 

Heatex 

Trise  =  8. 160x10  1.080x10  1.232x10  ^V-t-8.570xl0-* 

Twisted  Tape 

Trise  =  4.070x10"*V^-»-4.451x10"'V2+1.711x10^V-6.440x10^ 

where:  Trise 
V  = 

=  temperature  rise  (K) 
fluid  velocity  (m/s) 
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Figure  B.3  Friction  Temperature  Rise  Curves  for  Heatex  Insert,  Wire  Wrap 
Insert,  Twisted  Tape  Insert,  and  No  Insert. 
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Flow  Velocity  (m/s) 


APPENDIX  C.  SYSTEM  INTEGRITY  /  LEAK  TESTING 

As  mentioned  in  section  III  C,  the  material  for  the  auxiliary  condenser 
penetration  plates  was  changed  from  aluminum  to  stainless  steel  on  24  January  1992. 
The  stainless  steel  screw  thread  connectors  in  the  aluminum  side  plate  had  loosened 
to  the  extent  that  leakage  could  be  detected.  The  cause  of  this  loosening  was  due 
to  thermal  cycling  of  the  apparatus  and  the  differential  contraction/expansion  of  the 
aluminum/stainless  steel  combination. 

An  initial  leak  test  was  conducted  from  20  December  1991  through  2  January 
1992;  the  results  are  shown  in  Figure  C.l.  The  initial  mean  leak  rate  was  3.4  mmHg 
per  day. 

Subsequent  to  the  structural  modification  noted  above,  another  leak  test  was 
conducted  6-19  February  1992;  the  results  are  shown  in  Figure  C.2.  The  mean  leak 
rate  was  found  to  be  1.7  mmHg  per  day,  an  noticeable  improvement.  In  general,  a 
leak  rate  of  2  mmHg  per  day  is  considered  acceptable. 
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Figure  C.l  Apparatus  Leak  Test  1 
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Figure  C.2  Apparatus  Leak  Test  II 
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APPENDIX  D.  COMPARISON  OF  nXED  C,  vs  FLOATING  C.  SOLUTION 
METHODS  FOR  MODinED  WILSON  PLOT  DATA  REPROCESSING. 

As  related  in  section  VI  D  the  "fibied"  C,  method  and  "floating"  C,  method  were 
evaluated  against  the  original  instrumented  tube  results  in  order  to  choose  the  most 
accurate  method  of  reprocessing  non-instruroented  data.  The  results  of  the  fixed 
coefficient  versus  modified  Wilson  plot  floating  coefficient  method  comparison  are 
shown  in  Table  D.  The  fixed  coefficient  method  was  determined  to  be  more  accurate 
than  the  floating  coefficient  method  75%  of  the  time.  For  the  three  coolant  flow 
rates  considered  for  each  run,  the  overall  fixed  coefficient  method  mean  error  was 
±2.0%,  and  the  floating  coefficient  method  mean  error  was  ±5.4%. 

For  Re  <20,000  the  error  was  noticeably  higher,  which  tends  to  support  the 
assertion  of  a  number  of  researchers  that  data  with  Reynolds  numbers  below  20,000 
should  not  be  used,  since  the  flow  may  not  be  jfully  turbulent. 
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Table  D.  COMPARISON  OF  FIXED  C,  vs  FLOATING  C,  REPROCESSING 
METHODS 


File  Name 

Flow 

Vel. 

Film 

AT 

Inst 

Tube 

Data 

Fixed 

Method 

% 

error 

Float 

Method 

% 

error 

v, 

(m/s) 

ho, 

(kW/m^K) 

ho 

ho 

FIMAVSHl 

4.35 

29.33 

12.13 

11.93 

-1.7 

13.19 

8.7 

2.21 

22.81 

13.18 

12.82 

-2.6 

15.17 

15.1 

1.16 

18.27 

13.68 

13.16 

-3.8 

17.20 

25.7 

-2.7* 

+  16.5 

FIMAVSH2 

4.34 

36.09 

11.10 

11.14 

0.4 

11.07 

0.3 

2.22 

30.39 

12.03 

11.87 

-1.3 

11.72 

-2.6 

1.16 

24.62 

12.66 

12.50 

-1.3 

12.26 

-3.1 

-0.7* 

-1.8 

FIMAVSH3 

4.35 

46.10 

10.31 

10.26 

-0.5 

10.62 

3.0 

2.23 

39.35 

10.89 

10.83 

-0.6 

11.44 

5.1 

1.16 

32.56 

11.42 

11.40 

-0.2 

12.43 

M. 

-0.4* 

+5.6 

FIMAVSH4 

4.35 

55.20 

9.263 

9.146 

-1.3 

9.584 

3.5 

2.23 

47.31 

9.869 

9.756 

-1.1 

10.52 

6.6 

1.16 

39.45 

10.41 

10.39 

-0.2 

11.71 

12.5 

-0.9* 

+7.5 

FIMAVSH5 

4.35 

32.83 

8.618 

8.478 

-1.6 

9.018 

4.6 

2.23 

28.05 

9.589 

9.341 

-2.6 

10.36 

8.0 

1.16 

22.99 

10.23 

9.865 

-3.6 

11.66 

14.0 

-2.6*  +8.9 
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File  Name 

Flow 

Vel. 

Film 

AT 

Inst 

Tube 

Data 

Fixed 

Method 

% 

error 

Float 

Method 

% 

error 

V* 

(m/s) 

(kW/m^K) 

bo 

ho 

FIMAVSH6 

4.35 

39.23 

9.073 

8.953 

-1.3 

9.084 

0.1 

2.23 

33.76 

9.757 

9.447 

-1.3 

9.667 

1.0 

1.16 

27.98 

10.18 

10.25 

0^ 

10.64 

15 

-0.7* 

+  1.9 

FIMAVSH7 

4.35 

49.04 

8.960 

8.939 

-0.2 

8.637 

-3.6 

2.22 

41.25 

9.685 

9.513 

-1.8 

9.010 

I 

0 

1.16 

34.45 

10.17 

10.01 

-1.6 

9.209 

-9.4 

-1.2* 

-6.7 

FIMASW3 

4.35 

55.08 

9.379 

9.384 

0.1 

9.223 

-1.7 

2.22 

45.50 

10.15 

10.11 

-0.4 

9.813 

-3.3 

1.16 

36.67 

10.73 

11.18 

12 

10.62 

-1.0 

+  1.3* 

-2.0 

FIMASW4 

4.35 

52.83 

9.531 

9.558 

0.3 

9.528 

0.0 

2.22 

44.76 

10.15 

10.22 

0.7 

10.17 

0.2 

1.16 

36.16 

10.67 

11.18 

18 

11.08 

12 

+  1.9 

+  1.5* 

FIMASW5 

4.35 

54.28 

9.405 

9.453 

0.5 

9.367 

0.4 

2.22 

45.39 

10.08 

10.14 

0.6 

9.985 

0.9 

1.16 

36.69 

10.61 

11.20 

16 

10.90 

17 

+2.2  +1.3* 
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File  Name 


FIMASN4 


FIMASN6 


Flow 

Film 

Inst 

Vel. 

AT 

Tube 

Data 

V. 

(m/s) 

(kW/m^K) 

4.35 

46.24 

9.947 

2.22 

34.27 

11.12 

1.16 

24.82 

12.21 

4.35 

45.93 

9.872 

2.22 

34.40 

10.84 

1.16 

25.52 

11.81 

Fixed 

% 

Float 

Method 

error 

Method 

ho 

ho 

10.07 

1.2 

9.512 

11.34 

2.0 

10.24 

-j-1.6* 

9.823 

0.5 

9.302 

10.86 

0.2 

9.898 

14.35 

21.5 

-*-7.4 

11.77 

% 

error 


-4.4 

-7.9 

-6.2 

5.8 

8.7 

03 

-1-4.9* 
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APPENDIX  E.  UNCERTAINTY  ANALYSIS 


In  the  measurement  of  a  physical  quantity,  there  will  always  be  a  difference 
between  the  measured  value  of  the  quantity  and  the  actual  value.  The  magnitude 
of  this  difference  depends  on  the  accuracy  of  the  measuring  device  calibration, 
operator  experience,  environmental  effects,  etc.  Eventhough  the  error  associated 
with  a  single  measurement  may  be  rather  smaU,  the  error  may  grow  to  substantial 
proportions  when  combined  with  other  measured  quantities  in  a  given  calculation 
scheme.  The  best  estimate  of  the  difference  between  a  calculated  or  measured 
quantity  and  the  actual  value  o^  the  quantity  is  known  as  the  uncertainty. 

The  uncertainty  may  be  estimated  by  the  method  of  Kline  and  McCUntock 
[Ref.  24].  This  states  that  for  a  quantity  R,  which  is  a  function  of  several  measured 
quantities  (R=R(x,,X2,X3,...,x„)),  the  uncertainty  in  R  is  given  by  the  following 
relation: 


IP. 


\2 


w. 


Sx,  ’ 


1/2 


(E.1) 


where: 

Wr  =  the  uncertainty  of  the  desired  dependent  variable 
x„X2,X3,...,x„  =  the  measured  independent  variables 
Wi,W2,W3,...,W„  =  the  uncertainties  in  the  measured  variables. 


Georgiadis  [Ref.  5]  gives  a  complete  description  of  the  uncertainty  analysis  for 
this  experiment.  The  uncertainty  analysis  program  written  by  Mitrou  [Ref.  25]  was 
used  to  calculate  uncertainties.  The  uncertainty  for  runs  using  an  insert  ranged  from 
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±2.2%  to  ±4.6%.  The  uncertainty  for  runs  without  an  insert  ranged  from  ±10.9% 
to  ±17.1%.  Sample  ou^uts  of  the  uncertainty  evaluations  are  included  in  this 
appendix. 
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DftTft  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  FNMA'.'SHl 

Preeeure  Condition:  Uacuun 

Uapor  Temperature 
Water  Flow  Rate  (X) 

Water  Uelocity 
Heat  Flux 

Tube-metal  thermal  conduc. 
Sieder — Tate  conatant 


S3.7E4 
50.00 
2 . 7E 

3.7E5E+05 

385.0 

0.2200 


<Deg  C) 
(m/a  ) 

(  ^  7  ) 

(  W  /  m .  K  ) 


UNCERTAINTY  ANALYSIS : 


I  ItkO  T  ADI  c 

VI  I  I  \  A  I  I  ^  ^ 


PERCENT  UNCERTAINTY 


R  w  ^  c  i  w  *  N  i,*  ^  P  * 

Hcst  q 

Lcq-Nccn-Tcr;  Diff  LNTD 
u‘2i;  Rciistsr^cs,  Ru; 
Ovsrsil  .  T  ^  Uc 
Uq-tsr-Sidc  H.T.C.  Hi 
Urqcr-Sidc  ^.T.C.  Ho 


1  oc 

•  •  ^  w 

!  .4E 
1 .42 
.49 
4.40 
!  .50 

1  _  53 

2 .  E0 


pr.p  Tuc 


MMr  5  C  T  A  T  M  T  V  AKI  Al_  y  C  T  C  ; 


CMMAljCU  1 

1  1  .  ^  • . .  •  ^ 

- 

c 

3.567 

/  r 

.  wy  w 

- 

D/7l 

w  w 

(7i;7i 

*  w  w 

A 

tn 

(  ,/  «  ) 

■i 

4 

1  OOC  +  PiC 

.  <  W  W  W  ^  w  >./ 

{ W/>!'  Z  ) 

M-  •  ^  w 

Vpiv^  w  . 

3S5 

.0 

- 

0 

“•■jfTirTi 

>  b.  A.  W  ^ 

UNCERTAINTY  ANAL  V  5  S  • 


MAC  T  :£^c-i_r 


PERCENT  UNCERTAINTY 


Wall  Rcaistance,  R 


1  1 


H.T.! 


Urn* 

I  I  m 

X  Sh 


Side 


S  i  w  w 


H.T.C. 

H.T.C. 


I  Mjn 


u . 

<  •  A 

Ho 


0 

1 

1 

4 

1 

1 

•7 


!4 

es 

40 


31 

23 
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PQP  Tuc  UNCERTAINTY  ANALYSIS; 


File  Name ;  FNMAVSH 1 

PressLire  Condition:  Vacuum 


Vapor  Temper 

*  ^  «  I  ^  A 

W  «•  Uli  C 

- 

52 

C  1  t 

/  Deg 

Water  Flow  R 

m  f  y  \ 

Q  b  W  K  t»  t 

- 

30. 

v  %/ 

Water  Veloci 

4>  X, 

•'  r 

- 

1 

70 

(  iT  /  *  ) 

Heat  Flux 

m 

3 . 

235E+05 

(  W/r;'’2  ) 

imljC  iiCbSl  V 

hernal  conduc. 

m 

385 . 

(W/r^.K  ) 

Siede: — Tate 

conatant 

- 

n 

w  * 

^  A.  w  «y 

UNCERTAINTY  ANALYSIS ^ 


VARIABLE  PERCENT  UNCERTAINTY 


M-  <  .  p  t  S'  ■  0-4  «  MW 

2 

.  05 

Reynold*  Number  Re 

2 

.  1  9 

Heat  Flux  q 

2 

.  15 

Log-Mean-Ten  Diff,  LMTD 

.  34 

Wall  Resietance  Ru 

4 

.40 

n,.--.  1  1  u  T  r  1 1„ 

“  •  ^ . .  •  . 

.15 

Water-Side  H.T.C. ,  Hi 

,  00 

Vapor-Side  H.T.C.  Ho 

.  4  1 

nATA  pnp  Tuc  IJMPCCTATMTy 

£iMA^YS  I  ^  r 

File  Name:  FN.MA'‘^H4 

101  f  Rr  ‘ 

T*. — 

-W..KW'WVW.'b. 

s  1  f>iQC 

•  V  V  >  V'  w  w 

r 

'  wwy  w 

Flew  Rste  ••  %  ' 

»  90 

V'  w 

Ur t  er  Uelcc : t  y 

*  4 

*4 

( '  5  ) 

U«  •  *  C  "l  . . , 

‘  1  W  ^  b  '  Ah*'. 

w  7 

259E+05 

(  W  /  .-/2  ) 

•  395 

a 

(  W  /  .  K  ) 

SiwCjw^'"T*^e 

m  Ot 

%/ 

79Pif^ 

^ 

uncertainty  analvs is j 

t  lAD  T  Api  C 

V  I  •  I  '  A  t  I  w  k. 


PERCENT  UNCERTAINTY 


Maes  Flou  Rate  ,  Md 

0  90 

Reynolds  Number,  Re 

I  .12 

Heat  Flux,  q 

1  .00 

Leg-Mean-Ten  Diff,  LMTD 

*50 

•  w 

Wall  Resistance,  Rw 

4.40 

Overall  H.T.C.  ,  Uo 

1  .08 

Water-Side  H.T.C.  ,  H: 

1  .33 

Vapor-Side  H.T.C.,  Ho 

2.2E 
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DftTft  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  FNMAUSH4 

Pressure  Condition:  Atmospheric  (101  kPa ) 

Vapor  Temperature  ”  99.778  (Deg  C) 

Water  Flow  Rate  ( X )  -  50.00 

Water  Velocity  -  2.75  (m/s) 

Heat  Flux  -  B.713E+05  (W/m^Z) 

Tube-metal  thermal  conduc.  -  385.0  (W/m.K) 

Sieder — Tate  constant  “  0.2200 


UNCERTAINTY  ANALYSIS: 

VARIABLE  PERCENT  UNCERTAINTY 


Mass  Flow  Rate ,  Md 
Reynolds  Number,  Re 
Heat  Flux  ,  q 
Log-Mean-Ten  Diff,  LMTD 
Wall  Resistance,  Rw 
Overall  H.T.C.  ,  Uo 
Water-Side  H.T.C.  ,  Hi 
Vapoi — Side  H.T.C.,  Ho 

DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


(Deg  C) 

/  |yi  ^/  ■  ^ 

(W/m^Z  ) 
(W/m.K  ) 


D  •  •  •  •  • 

t  I  W  W  W  Ul  • 


Temperaturs 


Water  Flow  Rat 
Water  Velocity 
Heat  Flux 


1  •  !  V  \ 

V  W  '1*1 


Heat  Flux 

Tube-metal  thermal  conduc. 
Sieder-Tate  constant 


20.00 

1.15 

5.70BEH 

385.0 

0.2200 


1  _  27 
1  .49 

1  in 

I  •  W  t 

•  I 

4.40 
1  .40 
1  .55 
2.64 


UNCERTAINTY  ANALYSIS: 

IIAOTAOI  C 

V  it  •  «  A  t  «  W  4. 

Mas:  Flow  Rate ,  Md 
Reynolds  Number,  Re 
Heat  Flux  ,  q 
Log-Mean-Tem  Diff,  LMTD 
Wall  Resistance,  Ru 
Overall  H.T.C.  ,  Uo 
Water — Side  H.T.C.,  Hi 
Vapoi — Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

3.00 
3.11 
3.04 
.13 
4.40 
3.04 
2. 68 
4.57 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name: 

FNMAUSHB 

Pressure  Ccriditicn. 

Uacuum 

Uapor  Temperature 

- 

7E.243 

( Deg  C ) 

Water  Flow  Rate  (%) 

m 

80.00 

Water  Uelccity 

- 

4.35 

( rn  /  *  ) 

Heat  Flux 

- 

5.053E+05 

( U/m''2  ) 

Tube-metal  thermal  c 

endue. 

385.0 

(W/m.K ) 

Siede; — Tate  constant 

- 

0.2200 

UNCERT^^INTY 

DAD  T  AQI  c 

Vi  l>  V  A  I  I  k.  U. 


ANftLYS ! S • 


Flcui  Rste  Md 

I'WjrllWJikJW  •OktllLvCi  ,  I'W 


PERCENT  UNCERTAINTY 

e.S0 


1  1  /> 
I  a  >  W 

1  PkO 

I  t  W  W 


Lc5~Ncan-Ter^  D 
Wall  Raaiatanc 
O'a'crali  H.T.C. 


ff  LNTD 


Uiat 

cr — Side 

u 

.  T 

r 

I  W  •  y 

Hi 

^.^-C  ^  Am. 

U 

T 

r 

vwpw.  WkUl—  :c...W., 


40 

?'5 


DATA 


FOR  THE  UNCERTAINTY  ANALYSIS 5 


File  Ne*^e* 

CMMAijC 

HE 

Cerd: 

VAW>>*  vkjw-kiwi- 

V  .w  «Kw. 

•  m 

- 

-7C  C/TiO 

1  W  •  W  W  w 

(  n 

V  ^ 

Urn  *■  mr 

e  ■'  %  ' 

* 

zr>  piOi 

W  w  •  wv 

i.i  a  4  A..  0  4  V. 

m 

2 . 7£ 

( n  / 

Heat  Flux 

- 

4 . 707E+05 

( W/m 

^^4naik~*wNai4al 

rma 1  conduc . 

- 

•7C3C  Pi 

«./  w 

( W/m 

Siede^”T*^e  co 

Detent 

- 

%/  a  a.  W  W 

UNCERTAINTY 


ANALYS I S • 


\  lAD  T  AOI  C 


M  «  a 


Nd 

Qm 


f  -'w 

C  1  .  .  V  ^ 

IIWUV  I  A  f  kf 

Lcg-Mear.-Tcm  Diff  ,  LMTD 

IiIaI  1  Di.f 

«Vfa,IAA  /\WW4«VVW<'*k>C)  >>HV 

Overall  H.T.C.  ,  Uo 
Water-Side  H.T.C. ,  Hi 

He 


1 1 


u  T  r 


PERCENT  UNCERTAINTY 

1  .2E 
1  .48 
1  .40 

■70 

•  w  k/ 

4.40 
1  .4S 
1  .54 
2.62 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS 


File  Name:  FNMAUSHB 

Pressure  Condition:  Uacuun 

Uapor  Temperature 
Water  Flow  Rate  (X) 

Water  Uelocity 
Heat  Flux 

Tube-netal  thermal  conduc. 
Siedei — Tate  constant 


7B.541 
20.00 
1  .IB 

3.9BBE+05 

385.0 

0.2200 


( Deg  C  ) 

( m/s  > 
(W/m''2  ) 
<W/m.K  ) 


UNCERTAINTY  ANALYSIS: 

UARIABLE  PERCENT  UNCERTAINTY 


Hass  Flow  Rate,  Md 

2.99 

Reynolds  Number,  Re 

3.09 

Heat  Flux  ,  q 

3.03 

Log-Mean-Tem  Diff,  LMTD 

. !  9 

Wall  Resistance,  Rw 

4.40 

Overall  H.T.C.  ,  Uo 

3.04 

Water-Side  H.T.C.  ,  Hi 

2. BE 

Uapor-Side  H.T.C.  ,  Ho 

4.55 

DATA 

FOR  THE  UNCERTAINTY  ANALYSIS: 

C  «  1  « 

>  A  A  Nw 

Name : 

CKJMACQC  1 

« 

•  I  W  W  W 

ure  Condition: 

Atmospherio  (101  kPa  ) 

1  1  a  m 

V  I_1  p  NW  < 

-  100.138 

(Deg  C) 

Water 

Cl  .  Da  +  a  {  V  \ 

1  A  'm/y»  I'kWVW  \/*/ 

”  80 . 00 

Wat  er 

IJa  1  +  w 

VWAWWA  ^  Jf 

“  4.32 

( m/s  ) 

Heat 

Flux 

-  1.291E+0E 

(W/m'2  ) 

Tube- 

metal  thermal  oo 

nduc.  •  385.0 

(W/m.K  > 

Sieder — Tate  constant 

-  0.2200 

UNCERTAINTY  ANALYSIS: 

UARIABLE  PERCENT  UNCERTAINTY 


Mass  Flow  Rate ,  Md 

0.81 

Reynolds  Number,  Re 

1  .IE 

Heat  Flux  ,  q 

oc 

•  ww 

Log-Mean-Tem  Diff,  LMTD 

•JO 

Wall  Resistance,  Rw 

4.40 

Overall  H.T.C. ,  Uo 

.98 

Watei — Side  H.T.C.,  Hi 

1  .3B 

Uapoi — Side  H.T.C.,  Ho 

2.30 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  FNMAF051 

Pressure  CorcJition:  Atmospherx 

c  (101  kPa) 

ww^wi  •  wi'i^wt  Qtwi  C 

100.116  (Deg 

C) 

Water  Flow  Rate  (X) 

50.00 

Water  Uelocity 

2.73  ( n/s  ) 

Hest  Flux  ■ 

1.133E+06  (W/m'Z 

) 

Tube-metal  thermal  ccnduc. 

395.0  (U/m.K 

) 

Siede: — Tate  constant  - 

0.2200 

UNCERTAINTY  ANALYSIS^ 

ij^p T Agi_c  PERCENT  UNCERTAINTY 

Mass  Flow  Rste .  Nd 

1  ^  2S 

Reynolds  Number ,  Re 

1 .54 

Hest  Flux  q 

!  .  36 

Lc3“N£ar“'^ffri  Diff  LNTD 

1  c 

.  >  w 

l.l*  1  1  D*«  1  4  O.  . 

wwjiA  ,  I'w 

4.40 

Overall  H . T . ^ “ o 

1  ^  37 

).)« <1.  u  T  r  u^ 

y.  w  .  wi  «  *  .  •  4. 

1  , 59 

U  T  r  Ur^ 

yWMW-  W4«-<N.  l.w 

2 .  E  9 

n^TA  pop  Tuc  ijf (PcpTArfoTv  ANALYSIS: 

r  ,  rMMACiTici 

-  '’01  1 P  a  1 

*v'"pc*~  T*rq£r*tu'"£ 

1  ,71*7  1  /  Ha^ 

>  V  \<  •  W  •  <  ’  w  w  ^ 

C  ) 

£  t  £  ■*  F  1  c  L.'  R  a  t  £  '  %  ■ 

r  w 

«.  C'  •  V  «./ 

l.l»  y  U«  1  ^  i  ..  M 

1  . 1  E  (  n  /  s  ) 

H*£^F1u^  *” 

9.513E+05  (W/n'’2 

) 

395.0  (W/m.K 

) 

r>r^«4am4  v 

0.2200 

MMTCDTATMTV  AMAl  VCTC* 

uApTAD^c  PERCENT  UNCERTAINTY 

M«  •  •  Cl  ^.  .  O  m-^  m  Mr4 

1  «WM)  •>WvV«,  I'VJ 

3.01 

O4B 

.^w»-.ww  ,  .'W 

■7  1  T 

w  •  1 

Heat  Flux  ,  q 

3.05 

1  AM.  n«  tf  1  MTn 

«C''  1  f 

.09 

U'all  Rssislarcs  Ru; 

4.40 

Overall  H.T.C.  ,  Uo 

3.05 

Water-Side  H.T.C. .  Hi 

2.70 

vu^wt  wamC  '••'•w*,  i'nJ 

4.59 

• 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name: 

Pressure  Condition: 
Uapor  Temperature 
Water  Flow  Rate  (X) 
Water  Uelocity 
Heat  Flux 

Tube-metal  thermal  c 
Sieder-Tate  constant 


FNMAF20t 

Atmospheric  (101  kPa > 

-  99.787 

-  80.00 

4.33 

-  1.431E+06 

onduc.  -  385.0 

0.2200 


(Deg  C) 

(m/s  ) 
(W/m''2  ) 
(W/m.K  > 


UNCERTAINTY  ANALYSIS: 


UARIABLE 


Mass 

1  A  wU  i^aX  w  1 

Md 

Reyn 

olds  Number 

.  Re 

Heat 

Flux,  q 

Log- 

Mean-Tem  Di 

« t  1 

II,  w. 

MID 

Wall 

Res i stance 

.  Ru 

Over 

all  H.T.C.  , 

Uo 

Wate 

r-Side  H.T. 

C.  .  H 

A 

Uapo 

r-Side  H.T. 

r  u 

W  «  ,  II 

c 

PERCENT  UNCERTAINTY 

0.80 
1.14 
.  95 
.20 
4.40 
.  97 
1  .34 
2 . 29 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS! 


File  Name:  FNMAF201 

*  '  W  WW>>W*k*W*«» 

Uapcr  Tempers t ure 
+  F'c'  p-*-  /  •/ ' 

W  *  V  w  p  V  *  1  c  c  i  V  y 

U*.4  Cl 

■  <  C  U  V  <  A  M 

Tube-metal  thermal  conduc. 
Sieder-Tate  constant 


.  -  /  1 W 1  W  D,  \ 

AW  V  I  V  >  IV t  >a  « 

100.075 

50.00 

2.74 

9.941E+05 

TOC  Ql 

W  W  ^  •  lU 

0.2200 


(Deg  C> 

( m/s  ) 
(W/m''2  > 
(W/m.K  > 


UNCERTAINTY  ANALYSIS: 

ilADTADI  r 

V  (l  I  \  A  I  I  w  W  W 


Mass  Flou!  Rate ,  Md 
Reynolds  Number,  Re 
Heat  Flux  ,  q 


Log-Mean-Tem  Diff, 

LMTD 

Wall  Resistance,  R 

Iw 

Overall  H.T.C. ,  Uo 

Water-Side  H.T.C.  , 

Hi 

Uapor-Side  H.T.C. , 

Ho 

PERCENT  UNCERTAINTY 

1  ■?’7 

•  •  • 

1  .52 

1.35 

.14 

2.B7 

1  .36 

1  .57 

2.66 
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DftTft  FOR  THE  UNCERTftlNTY  ANALYSIS 


File  Nane: 

Pressure  Condition: 
'v'apcr  Temperature 
Water  Flow  Rate  (X) 
Water  Velocity 
Heat  Flux 

Tube-retal  thermal  c 
Siede; — Tate  constant 


FNMAF201 

Atmospheric  (101  kPa  ) 

-  100.142 

-  20.00 

-  1  .  IS 

-  E.405E+05 
nduc.  "•  395.0 

0.2200 


(Deg  C> 

( m/s  > 
<W/m'2  ) 
(W/m.K  > 


UNCERTAINTY  ANALYSIS: 

nApTjagi_£ 


PERCENT  UNCERTAINTY 


mass  I- 


D.+ 


^ate  ,  Md 

.  .  w  w-  t.  I  *  w  H 


Ptff  I  MTn 


!  er.  Uit 

i.l»  11  Oa.  .  «  4  A  o. 


n  w .  -  1  1  u  T  r  u  ^ 
L1.4  ••--C  ^  Am  u  T  r 

V*w  W.  -  *  .  •  .  .  .  w  .  , 

1  y  Am  u  r  r 


3.01 


1  T 
>  ^ 

05 

09 

G7 

OkC 

w  ^ 

nrh 

•  ’L/ 

CQ 


FOR  juR  UNCERTAINTY  ANALYSIS t 


C  ,  \  m  f^\mm,m  • 
t  A  «  W  I  .  w  I  I  W  • 

v'Spor 

Z  ]  .  Cm*  m  .  ®/  ' 

<NwVW'  AWW  ‘'WVW 

m-  ilm'i  ^ 

y.-.vw 

Umm*  Cl.., 


FSONNASN ] 


T.  _  _ 

.  >.«  C.I  C  I  '  V  A 


,  1 


3UC  . 


•  <  1  0  1  ) 
i  OtOi  t  cc 


549E'i-05 

0 

0130 


i  Dm 


{  W/rs 
(W/m 


r  I 


»  ) 

2  ) 


UNCERTAINTY  ANALYSIS: 

HAD  T  AQ|_C 

PERCENT  UNCERTAINTY 

Mass  Flew  Rate,  Md 

0.S0 

Reynolds  Number,  Re 

1.11 

Heat  Flux,  q 

1  OiO 

i  •  v* i. 

Log-Nean-Tem  Diff,  INTO 

.45 

Wall  Resistance,  Rw 

2.E7 

Overall  H.T.C.  ,  Uo 

1.11 

Water-Side  H.T.C.  ,  Hi 

15.41 

Uapor-Side  H.T.C.  ,  Ho 

10.91 
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DftTft  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  FSONMASNl 

Pressure  Condition:  Atmospheric  (101  kPa  ) 


Vapor  Temperature 

- 

99.842 

( Deg 

Water  Flow  Rate  (X) 

- 

50.00 

Water  Velocity 

- 

■?  ■7C 

^  •  1 

( m/s  ) 

Heat  Flux 

- 

4.012E+05 

(W/m-Z  ) 

Tube-metal  thermal  conduc. 

- 

395.0 

(W/m.K  ) 

Siedef — Tate  constant 

- 

0.0130 

UNCERTAINTY  ANALYSIS: 

VARIABLE 

PERCENT  UNCERTAINTY 

Mass  Flow  Rate ,  Md 

1  _  27 

Reynolds  Number.  Re 

1  .50 

Heat  Flux,  q 

1.39 

Log-Mean-Tem  Diff,  LMTD 

.32 

Wall  Resistance.  Rw 

0  C-7 

^  •  W  1 

1  1  u  T  r  11,. 

1  .42 

Water-Side  H.T.C.  ,  Ki 

15.44 

1  u  T  r 

WAWW  ri«l*W*f  IIW 

17.07 

na|T|a  cgp  juc 


UNCERTAINTY  ANALYSIS: 


C  ^  \  m.  K\m^m  * 

FSONMASNl 

<  WWWW<  W  w'W>iW*VAw.<* 

Atmospher 1 

c  (101  hPa ) 

iwiapC>  WWW  w 

- 

1 p»7C 

WWIIrW<  1  AWW 

- 

70.00 

WUVWI  VWAWWAV/ 

■i 

3.81 

Heat  Flux 

- 

4.39BE+05 

Tube-metal  t herma 1  c 

onduc .  “ 

395.0 

Sieder  Tate  constant 

- 

0.0130 

(Des  C) 

(  m  /  s  > 

(  y  ,/(*.'■  2  ) 
(W/m.K  ) 


UNCERTAINTY  ANALYSIS: 


I  lAD  T  AP<  C 

V  t  iiv  A  n w w 


PERCENT  UNCERTAINTY 


Mass  Flow  Rate ,  Md 

0.92 

Reynolds  Number.  Re 

1  .21 

Heat  Flux  ,  q 

1  .09 

Log-Mean-Tem  Diff,  LMTD 

.41 

Wall  Resistance,  Rw 

2.67 

Overall  H.T.C.  ,  Uo 

1.17 

Watei — Side  H.T.C.  ,  Hi 

15.42 

Vapoi — Side  H.T.C.,  Ho 

12.42 
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APPENDIX  F.  INSTRUMENTED  TUBE  CONSTRUCTION 


The  medium  diameter  instrumented  tube  of  Poole  [Ref.  6],  with  six  waD 
thermocouples,  was  fabricated  by  welding  together  a  copper  tube  in  three  pieces 
after  imbedding  capillary  tubes  at  mid  depth  in  the  tube  wall.  Teflon  thermocouples 
were  inserted  into  the  capillary  tubes  for  tube  wall  temperature  measurement. 

An  attempt  was  made  during  this  study  to  construct  instrumented  tubes  of 
small,  medium,  and  large  diameter  in  the  following  steps: 

1.  Take  thick-waUed  copper  tube  bar-stock  and  machine  to  specified  inside 
diameter,  outside  diameter,  and  length. 

2.  Cur  four  evenly  spaced  slots,  of  given  depth  and  width  to  accommodate  metal 
sheathed  thermocouples,  from  halfway  along  the  tube  to  the  end. 

3.  Solder  metal  sheathed  thermocouples  into  groves. 

4.  Cur  copper  strips  from  another  tube  to  fit  into  the  top  of  the  slot. 

5.  Clamp  copper  strips  to  the  slots  (using  jubilee  clips)  and  solder  in  place. 

6.  Turn  off  excess  copper  from  strips  to  original  outside  diameter. 

7.  Send  tubes  to  plating  shop  and  plate  the  whole  tube  to  a  given  thickness  of 
plate. 

8.  Return  tube  to  NFS  machine  shop  and  machine  back  down  to  desired  outside 
diameter. 

The  process  was  completed  through  step  2  of  the  above  procedure  and  the 
tubes  are  ready  for  step  3.  Several  of  the  latter  steps  were  attempted  using  a 
practice  tube  with  monel  wire  placed  in  the  slots  on  the  tube  vice  thermocouples;  this 
was  to  ensure  the  process  was  safe,  and  to  prevent  destruction  of  assets  (small, 
medium,  and  large  diameter  tubes  machined  to  specifications,  and  the  associated 


metal  sheathed  thermocouples).  The  following  summary  documents  the  steps  taken 
and  lessons  learned. 


1.  Completed  tube  fabrication  through  step  2  of  the  fabrication  procedure 
(above).  A  practice  tube  was  used  from  this  point  to  continue  the  procedure 

2.  Cutting  the  copper  strips  proved  difficult,  so  the  monel  (similar  to  the 
stainless  steel  sheath  on  the  thermocouples)  wires  were  silver  soldered  after 
being  pinned  in  place.  No  lead  was  permitted  in  the  solder  since  copper  will 
not  plate  to  lead,  but  will  plate  to  silver. 

3.  The  practice  tube  was  sent  to  a  local  contractor  for  electroplating.  (Bay 
Custom  Chrome,  Marina,  Ca.) 

4.  Two  different  plating  procedures  were  used  on  the  tube; 

a.  First,  a  cyanide  bath  treatment  was  used  to  electroplate  the  tube. 
However,  this  procedure  severely  scorched  the  surface  of  the  tube, 
resulting  in  the  return  of  the  tube  to  the  NFS  machine  shop  for  repair. 

b.  After  repair,  the  tube  was  again  sent  off  and  subsequently  treated  with 
an  acid  bath  procedure,  which  resulted  in  an  acceptable  tube  surface 
finish. 

5.  The  original  thickness  of  electroplate  was  not  thick  enough  for  our 
application,  so  another  acid  bath  treatment  was  performed  to  bring  the 
surface  thickness  to  the  desired  level. 

6.  The  final  product  was  suitable  for  our  intended  application  with  the  exception 
of  the  following  items  which  need  to  be  corrected: 

a.  Existing  voids  (from  non-uniform  soldering)  or  flaws  in  the  tube  surface 
would  not  fill  in  to  give  a  uniform  thickness  around  the  tube. 

b.  The  practice  wires,  tube  end  pieces,  and  tube  interior  were  not 
protected  properly  during  the  electroplating  process,  resulting  in 
unacceptable  copper  deposition  on  critical  areas. 


It  is  recommended  that  the  instrumented  tube  efforts  be  continued  with 
particular  attention  to  the  following  items: 


1.  Silver  soldering  is  an  acceptable  method  of  fixing  the  thermocouples  in  place, 
since  the  thermocouple  melting  temperature  is  ~1700°C  (silver  solder  is 
applied  at  ~1100°C). 
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2.  If  the  voids  left  in  previous  solderine  efforts  persist,  then  cut  copper  strips  to 
place  over  the  thermocouples  to  hold  them  m  place. 

3.  Ensure  that  critical  areas  on  the  tube  are  well  protected  during  the 
electroplating  process. 
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APPENDIX  G.  DRPINST  PROGRAM  LISTING 

The  program  DRPINST,  which  was  used  to  collect  all  instrumented  data,  is 
listed  in  this  appendix. 
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DRPINST 

WRITTFN  FOR  INSTRUMENTED  TUBES 
BY  S.  MEMORY  (i0TH  DECEMBER  !S9!) 

USED  TO  COLLECT  DftTft  ON  THE  INSTRUMENTED 
TUBE  FABRICATED  BY  POOLE  (1983)  FROM  WHICH 
NEW  INSIDE  HEAT-TRANFER  CORRELATIONS  WERE 
EMPIRICALLY  DETERMINED  FOR  NO  INSERT.  WIRE 

UDAp  TMCCCT  JMCCpT  CftCCC  gy 

u;  _  CyCMCCM  l  peg  t  QQ2  )  , 

All  instrument  calibrations  for  flow  METER 
TuroMQrQijpycc  THERMOPILE  PRESSURE  TRANSDUCER 
friction  ’’’EMP  rise  for  no  insert,  WIRE  WRAP  INSERT, 
HEATEX  INSERT  AND  TWISTED  TAPE  INSERT  WERE 
pccpnpMcn  gv  cycKjccKj  ( SEP~DEC  1991  ). 
ijADjnijc  ppngcAMC  yepe  yeen  tq  DgppQrecc  hat^ 

COLLECTED;  THE  APPLICABLE  FILES  WITH  THE  APPROPIATE 
c  gDpnr  r  c  c  T  Kjg  ppncp^MC  Ape  ac  FOLLOWS; 


0T  7  1 

files 

REPROCESS 

TKiC  DDHCDAM 

«  .'WWIvt... 

FLFNFFN  FTMA‘JFHI-7 

nppTMCT2 

.71  7  C  1 

CupMCCKj  CTMACM4  C 

npCTMCTI 

'?  *• 

pycyceM  e  j  MA._icy?_p 

npp  T  MCT  1 

,71"?^. 

r  i-‘  F  S  E  -  -  I  ^  F  ‘  p  y*.! 

T  c  r  Pi  A  T  A 

nppf.;  c 

,  A.  ,  c^:  C  T  1  C  C 

f?  “  4  . 

S  L'  T  T  E  N  D  0  F  ^  D  F  T  F 

ORPK  S2 

nTC  i 

,  JAM  OpTTCKj  flATA 

nppuco 

ny.c  . 

W  W 

j777  i 

FILF  DEFINITIONS 

cQP  SWENSpN 

: 

1 

04  01 

LE’^'IER  /NU^SER 

definition 

POS I T I  ON 

(  CT 


CONDENSATION  CONDITION 
D-DROPWISE,  F-FILMWISE 
TUBE  TYPE  I-INSTRUMENTED 

MmMriM-  T  MCTpl  (MCMTCn 

TUBE  TYPE  M-MEDIUM 
PRESSURE  CONDITION 
A-ATMOPHERIC  (101  kPA  ) 
U-UACUUM 

WHEN  A  &  U  APPEAR  TOGETHER  AS  IN 
FIMAUSH1  IT  REFERS  TO  A  SERIES 
OF  RUNS  AT  VARIOUS  PRESSURES  AND 
UAPOR  UELOCITIES. 

TUBE  TYPE  S-SMOOTH 
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1057! 

F-FINNED 

10581 

BTH  TYPE  OF  INSERT  INSTALLED 

1059! 

N-NONE 

10E0! 

U-UIRE  URAP 

1051  ! 

H-HEATEX 

1082! 

THE  HEATEX  INSERT  WAS  USED  EXCLUSIUELY 

10B3' 

FOR  FINNED  TUBE  DATA.  SO  POSITIONS 

10S4 ! 

B  &  7  REFER  TO  THE  TUBE  FIN  SPACING 

!0E5‘ 

05“0 . 5rnn 

106B  ! 

1 0“1 . 0npi 

10B7' 

15”1 .5nn 

1  0B9 ' 

20"2 . 0nn 

1  OiC  □  1 

1  W  W  u 

FINAL  POSITION  RUN  NUMBER 

1  070  1 

t  an  1  1 

1  Vf  •  < 

MEANING  OF  ALL  FLAGS  IN  PROGRAM 

1072  ‘ 

1  073  i 

IFT:  FLUID  TYPE 

1  074. 1 

ISO:  OPTION  WITHIN  PROGRAM 

1  anc  1 

•  b  <  ^ 

IM:  INPUT  MODE 

1  anc.  1 

>  «./  1  w 

TC0;  FINNED  OR  SMOOTH 

1  077  1 

INN:  INSERT  TYPE 

1078* 

lUT:  LOOP  NO.  WITHIN  PROGRAM 

1079* 

ITDS:  TUBE  DIAMETER 

1  oioa  1 

<  ^  w  w 

TOP-  POCCCIIPC  PPlMnTTTnKJ 

1  1  1 
•  u  w  • 

Ton*  niiTDHT  DcniiTocn 

WWII  W< 

1082 

COM  /Cc/  C(7) 

1  003 

COM  ^Cc5E/  T££(E) 

1  aoA 

•  V  W-  ^ 

COM  /Cc5£/  T58  (  £  !’ 

1  00C 

COM  'Cc£7/  T£7^  £  ) 

1  0QC 

COM  ■'CcES/  T£S'«£) 

1  007 

rnw  /F  1  ^  *  f  t 

1  099 

HTM  f  1  n.  \  T../CN  T.  /  C  \  T..WC\ 

^  A  ‘  1  ^ili> 

1  099 

QATA  0.10095091  2£727. 94389  -78734E.9295  792 

:)ccqc  01 

A.  ^  ^  W  •  W  1 

1  aQPi 
•  «/ 

HATA  -a'?A'740CCaQ  C  QTCOOCu.ll  _*?  CClQ7Cj.l'7  T 

9407eE+!4 

1  09 1 

READ  C(  ? 

1  ^^0“? 

1  k 

HAT^  273. 1£  2.£945E~2  “7.2S71E-7  3.2941E“11, 

-9.7719E-! 

1  aQ-x 

•  w  w  ^ 

READ  T5£(  • ) 

1  094 

DATA  273. 15,2. E879E-2 .-£. 88S3E-7 ,-3. 1242E-1 1 

,1 .3275E-1 

1  09E 

DCAn  TCC  (  *  ^ 

iV  W  •*.'  ■  W  '  '  ! 

1098 

DATA  273. !5,2.5823E-2  ,-7 . 3933E-7 ,2 . 8S25E-! 1  . 

1 .9717E-15 

1  aon 

1  ^  1 

pcan  TC7/  « ) 

!  099 

DATA  273.15.2.5931E-2  .-7 . 5232E-7 .4 .05B7E- 1 1  . 

-1 .2791E-! 

1099 

READ  T59' •  ) 

1  t  PlOl 

1  >  %y  %/ 

Dssp-.!524  '  Inside  dianeter  of  stainless 

steel  test 

1  101 

Ax-PI*Dssp''2/4 

1102 

L“. 13335  '  Condensing  length 

1  1  ffl-J! 

•  •  w 

L1-.0B0325  1  Inlet  end  “fin  length" 

1  104 

L2-. 034925  •  Outlet  end  "fin  length" 

1  1 

•  >  w  ^ 

PRINTER  IS  1 

1  10B 

BEEP 

6.9.7121E-20 

4, -1 .0198E-1S 
.-2.249BE-19 

5, G.4402E-20 
section 


109 


1  ffiQ 
>  ^  ^ 

1  1  0 
1  '  1 
1  1  2 
1  13 
1  14 
1  1  7 
1201 
1  T*  T 
1  2£ 

;  29 

1  y. 

t  7C 

1  3S 
1  4 1 

1  A  A 

]  47 
1  C0 

1  53 

1  cc 

'  7  ^ 

^  "7  / 


INPUT  " SELECT 
PRINT  USING  “ 
PRINT  USING  " 
PRINT  USING  " 
PRINT  USING  *' 
INPUT  Isc 


FLUID  (0-STEftM.  1-R-113,  2-EG)*‘,Ift 
4X/"*Select  opticn:“***‘ 

1  Take  data  or  re-process  previous  data*'"' 

2  Print  raui  data  (NOT  COMPLETE) . * 

3  Purge . 


cv  •• 

si;- 

BX 


A  3  <• 


THEN  30S2 


PRINTER  IS 

T  . 

A  J  W  «-»  V' 

CLEAR  709 

DCCD 


701 


TMpijT  "rfsiTcp  NONTH  date  AND  TIME  (  MMjDD'HHjMM;  SS  )"  Date$ 
OUTPUT  709i "TD" ;Date$ 


TC  T , 

4.  I  4  J  W  W 


1  THEN 

BEEP 

TMDMT  “CVTD  AMH  UTT  CMTCp*'  n|, 

CKJn  TC 


ENTER  709;DateS 

pCTMT  " 

PR  I  NT 

CC  T  MT  1  IC  7  Kit”  "  1  /7kV 

■  r  w>-.4'Mw  <  ■  y 


'  MHTC  • 


ate  and  tire  :";Date$ 

•  DRPINST* 


IF  I  ' 

DCCO 


^k«^  TUCK; 


INP''T  "ENTER  INPUT  MODE  «  1  "NEU  DATA  2*EXISTING  FILE)"  Ir 

DCCD 


1  c  a 


ruCM 


;  C  ' 

SEER 

1.(1 

tn 

I  f'i" '  'T 

'■CTnc  A  ^;AM£  rnc  p 

;  PE 

C  C  T  r.i  T 

jcTMC  ‘■■pile  "a'^e 

zo* 

r  EC ATC 

EDr"''  0  * :  1  e$  ■  ^ 

204 

ACC  T0KI 

@Rile  TO  0  files 

?(7i- 

*>  w 

BEEP 

1  .Ti 

INPUT 

•’CMTcc  iJEOME'^’P*  COO^ 

2  ‘  1 

PR'^NTER  IR  1 

?'  T  T 

I  rn-0 

T  1  C 

BEEP 

*?  1  Q 

PRINT 

"SELECT  INSERT  TYPE:" 

T  (71 

PRINT 

"  0  NONE  '^DEPAULT)*' 

2  2  1 

PR  I  NT 

"  1  "^UI^TED  TAPE" 

T  "7  7 

p  p  T  M  T 

'  7  I.ITDC  MCAD" 

to  W4l^^  mit.if 

"'•7"? 

PRINT 

-■  3  UCAJCV" 

225 

T  MDI  \T 

4  ■  W  1 

I  nn 

225 

OUTPUT 

«r . 1  - . T  Tp^ 

>  44Cl4»y,4r*l> 

7  7  "7 

T  MDI  IT 

4  •  W  ■ 

"NO.  OF  THERMOCOUPLES 

7  7  0 

T  C  1  fn 

4 .  4  .  g 

»1  TU£M 

229 

I NPUT 

"ENTER  FIN  HEIGHT  FIN 

*^30 

ELSE 

23 1 

Ck./Ti 

1  •  V' 

F!lE‘  D  ^lie$ 
:  ■  1  4  r  :  D  f  1 1  •  S 


1 -FINNED  0-PLAIN)"  If 


IN  WALL ".I nut 
PITCH  AND  FIN  WIDTH 


g 


(Fh.Fp.Fu) 


IN  MM“  .Fh  ,Fp  ,F 


no 


in  in  in  in  in  in  4:>  ^  4:^  4u  45.  4=«  i>j  i^i  i>i  i>i  i>j  im  i>j  i>i  ixi  w  ro  ro  ro  ro  fo  ro  ro  ro  f'J  to  •-  to  to  to  to  to  to  to  to  to 

4i  4:w  4^  —  (9  o  to  (D  CO  -o  -o  -o  -o  -o  in  in  in  in  in  in  -o  -o  in  in  in  in  m  in  in  4^  to  to  co  -o  -o  -o  cn  in  in  m  t"  in  in  4:^  4^  i>t  i>4  i>4  i>4 

ID  in  10  o  4^  —  CO  in  o  ID  to  -o  4:»  to  -0  cn  in  lo  to  CD  to  in  m  tn  i>i  o  -o  4:»  —  to  in  in  ro  cn  to  o  -0  4^  —  to  ti  in  to  10  cn  to  -o  in  4^  to 


Fp-0 
Fu,-0 
END  IF 

OUTPUT  eFilesInwt  .Fp.Fui.Fh 
ELSE 

IF  Ijob-1  THEN  1255 
BEEP 

INPUT  "GIUE  THE  NAME  OF  THE  EXISTING  DATA  FILE"  .D_f i le$ 

PRINT  USING  "1GX,""This  analysis  uas  performed  for  data  in  file 


,10A“;D_ 


IF  Ijob-1  THEN  1270 

Nrun- 1 8 

BEEP 

INPUT  "ENTER  NUMBER  OF  DATA  SETS  ( DEF-1 8 )"  .Nrun 

ASSIGN  ©File  TO  D_fileS 

ENTER  @File;IfQ.Inn 

ENTER  @File;Inut  .Fp.Fu.Fh 

END  IF 

IF  Ijob-1  THEN  1543 
PRINTER  IS  1 
BEEP 

ppjMT  "SELECT  TUBE  DIA  TYPE:" 


PRINT  "1  SMALL" 

poTMT  "2  MEDIUM  (DEFAULT)" 

PRINT  "3  LARGE" 

TMOIJJ  J  +  0« 

PRINTER  IS  701 

•i./ibj  4.W  w<  r^wUAUrt  W>>U  ItWwC* 

Do-. 01905  '  OD  of  medium  tube 


Do-. 01905  '  OD  of  medium  tube 

D1“. 01905  '  OD  of  unhested  inlet  length 

D2-. 015875  '  OD  of  unheated  cutlet  length 

Dr-. 015875  1  Thermocouple  burial  depth 


Di-. 009525 
Do-. 0127 
Dr-. 01 

CKin  TC 

Wl  ow  A  t 

IF  Itds-3  THEN 


ro  BE  MODIFIED  WHEN  KNOWN 


Dr-. 02 


!T0  BE  MODIFIED  WHEN  KNOWN 


Kcu-305 

Rm-Do*L0G( Do/Di  )/( 2*Kcu )  !  Wall  resistance  based  on  outside  area 
BEEP 

INPUT  "ENTER  PRESSURE  CONDITION  ( 0-V  .  1 -A  )"  . Ipc 
BEEP 

PRINT  USING  '16X,”"This  analysis  includes  end-fin  effect”"" 

PRINT  USING  " 1GX  “Thermal  conductivity  -  "",3D.D.""  (W/m.K) 

PRINT  USING  "1BX  .""Inside  diameter,  Di  -  "".DD.DD,""  (mm)"" 


""" jKcu 
" !Di*1000 


111 


1  c:c7 

•  W  4* 

PRINT  USING  "  IBX  " 

Oi^t 

1B00 

IF  Ijob-1  THEN  1648 

1S09 

BEEP 

\  C7C 

«  w  «/  w 

INPUT  "GIUE  A  NAME 

FOR 

1  039 

CREATE  BOAT  Wtf$.5 

1B42 

ASSIGN  @F:lsl  TO  Wtf$ 

’B48 

IF  Ijob-1  THEN 

t  CC  f 

T 

4  W  V  ^ 

1  BE4 

GOTO  1 702 

1  CC7 

>  w  «.>  • 

END  IF 

1  CCPi 

1  W  W  V 

BEEP 

1  093 

TMDMT  "SELECT  OUTPUT  ( 

1  CC  Q 

>  w  u  w 

TC  T  TUCM 

*•  ^  I'.i-i'i 

1  972 

PRINT  USING  "lEX 

Tub 

1  970 

ELSE 

1  979 

PRINT  U^IN‘“  "IBX  " " 

Tub 

1001 

PRINT  USING  " 1  Ex 

Fir. 

2 .  DD " 

.  c«  c. .  CU 

y  1  p  ,  •  w  ,  1  II 

7  994 

END  IF 

1  702 

T«i7i 

W  <4.‘ 

1  77  1 

/-vrs-1 

-  ,DD.DD. 


“  (  pf*m  V  " 


.dth,  or.d  Height  (nr*.):  "‘',DD.DD,2X 


7  nn  7v 

4.  ■  WW  ,  t 


774  c*p-2 1  >  ' 

777  7 • T i 1 


W  V' 

T  C  ' 

A  All 

1  TH^N 

793 

BEEP 

7S4 

nijTPU'^ 

7(70.*'  AD 

•  V  w  1  11' 

ftF40  4L40  UPE" 

-JuC 

OUTPUT 

70 Q .  *  AC 

C  A  " 

707 

I  nfut 

"  1  I  u  £  Tn 

I'uC."..  r.Kl  r  p.r*r-CNjT  O' AT 

7  C  C 

I  C  J  «  1 

T  u  c  f  j 

7G2 

0  U  T  ^  U  T 

700  .  "AC 

AC40  ‘^L41  “pP" 

-  QC 

0 UT^UT 

7i7lQ  .  ■  A  1 

w  ‘  « 

C  A  •’ 

73C 

END  IF 

■Da  T 

BEEP 

S04 

I N  P  U  T 

"CKiTCp  pi_ 

OUNETEP  PEf^DING"  Fn 

C07 

OUTPUT 

7(70  .  ■'  AD 

.  V  «  III' 

AC90  ALS2  UPS" 

c  *.  a 

W  i  V 

iTDl  JT 

w  w  1  1  W  1 

700  .  "AC 

SA" 

9  ■  3 

ENTER 

?09;Etp 

c  ^  9 

r\\  \TD\ \T 

W  W  1  I  W  1 

709 ;  *'  FiS 

C  A" 

C  1  u 

BEEP 

922 

I  NP‘  ‘T 

"  CONNECT 

U0LT4PE  LINE”  01^ 

C7C 

ENTER 

7  0  9 ;  E  V  c  1 

PI'S 

C0EP 

1 

I NPU*^ 

"DISCONNECT  UCLT.^GE  LINE'*  OK 

934 

TC  Ds,^ 

A  '  *  W 

!  ^  t  Turw 

940 

BEEP 

943 

I N  P  •  .*  T 

"  I NUAL I D 

nni  TAOC  -  Tpv  AOATM 

VWU.IilWk-  l~'i  i1wi1a«« 

946 

GOTO  1 

0  1  o 

w  >  ^ 

94  9 

END  IF 

^52  ' 

OMTpMT 

7(70  .  ••  AC 

1  w  1  1  1  w 

CA  " 

009  1 

CMTCD 

I « •  1 « 

709iUtrar 
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2002 
2009 
2008 
201  1 
2012 
2014' 
2015 
2017 


OUTPUT  709s “ftS  SA" 

ENTER  709sBaMp 
Etp-Etp*1 .E+B 

OUTPUT  709; "AR  AF40  AL53  UR5" 

Nn“7+Inwt 

FOR  1-0  TO  Nn 

OUTPUT  709; "AS  SA" 

IF  I >7  THEN 
Se-0 

FOR  K-1  TO  20 
ENTER  709 ;E 
Sc-Se+E 
NEXT  K 

Enf ( I  )-AB£(£e/20  ) 

ELSE 

ENTER  709 ;E 

!  T  \»ADC  t  C  \ 

I  N  A  >  iiww  V  w  / 

ENO  IF 

Er!f(  I  )-Epif  (  I  >*1  .E+6 
NEXT  I 


OUTPUT  709; "AS  SA" 
OUTPUT  713;"T1R2E" 

UA  T  T  ■> 

W  •  I  A  I 

ENTER  7 1 3 ; T 1 1 
OUTPUT  713;“ T2R2E “ 
WAIT  2 

ENTER  713;T2 
OUTPUT  1  7  ;  " 1  R2E  " 


ENTER  7 1 3 ■ T 1 2 

Tl«{Tl  \ mT ]  2  )*,£ 

OUTPUT  7 1 “  T3R2E  “ 


INPUT  “ENTER  PRESSURE  GAUGE  READING  (Pga)”,Pga 

Pvap1-Pga^5894.7  !PSI  TO  Pa 

ELSE 

ENTER  @Fi is; Bvcl  ,Banp  .Ptran  ,Etp,Enf(^),Fm,Tl ,T2  ,Pvap1 
IF  J-1  OR  J-20  OR  J-Nrun  THEN 


Ng-0 

TNin  TC 

^  I  At 

END  IF 

T5tear,1-FNTvsv57(Er,f  (0  )  ) 
Tstaani -Tsteanl -273.15 
T5tear*!2-FNTvsv5E(Enf<  1  >) 
Tst=an2-Tstean2-273. 15 
T5tean"(  Tstear.l  +T5teaP!2  )/2  . 
Tstean-Tsteami 
Troon-FNTvsv58( Enf ( 2  > ) 
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2018  Tcon-FNTvsv58(Emf(7)> 

2019  Troom-Troon-273 . 1 5 

2020  Tcon-Tcon-273. 15 

2021  Tu;Pi-0. 

2022  FOR  1-0  TO  Inwt-1 

2023  Twf I  )-FNTv5v57<Emf( 1+9  )  ) 

2024  Tij(  I  >-Tw(  I  )-273. 15 

2025  Tuim-Tun+Tw^  I  ) 

2029  NEXT  I 

2032  Twn-Tuw/Ir.wt 

2044  F 33 1 "FNP v5 1  (  Ta  t ear. ) 

2057  OUTPUT  709; “AR  ftF64  ALS4  UR5” 

2059  OUTPUT  709;  "AS  SA"  !  PRESSURE  TR.ANSDUCER 

20CQ  Ss»0 


oracpi 

V  w  W 


22B7 

1 

^  V  >  ‘ 

2072  ‘ 

2073 

9<^’7C 

4.  U/  I  W 

r»<7»'7'7 

«.  V/  •  • 


FOR  K-1  TO  20 

CMTCD  '7r^Q,C4^«r, 

Cii 

Ss-Sa+Etran 
NEXT  K 

s- /20  ) 

SEEP 

PRESSURE  IN  Pa  FROM  TRANSDUCER 
Pv3p2-(  -2 . 93G04^Ptrar,+  M  .7827  :‘*B994 . 7 
Pt eat 1-Pvap1  '  GAUGE 

p+e5+2-Pv3p2  ‘  TRANSDUCER 

CORRECTION  FOR  AIR  CONTAMINATION 
P^p-Pvap2^1 .E-3  f  TRANSDUCER  IN  kPa 
PI*  5*Paat •  1  , E“3  '  SAT.  PRESSURE  IN  kP* 


PKc 

T  « 


sOv.  A''  ^ 

I  •  V  O  t-i  I 

,  4  .CMT  .. 


!  .£-3 

D«  »4  \ 


J  CAMpr  TM  i:Vi 


2093 

2094 

'ifflQQ 

.•WWW 

7  t  J  /a 

4.  >  I  W 

2  M  3 
2  1  1  B 

7  1  1  7 
'M  1  Q 

W  :  *  w 

Z  1  34 
2 1  35 
2  1  37 

^  A 

*.  •  T  W 

*7143 
2  t  4E 
2149 
2  1  50 
]  52  ^ 
'y  1  cc 


T  . .  1  »CMTs.  •  r.  /  D  w  ».-  1  ^ 

.  «  c  .  II  I  .  WfW  •  •  V  wp  I 

T  , 'll  «  C  Ki  T  M  ^ 

.  r  LJK*-  •-<»■»  W^  •  ■  .  ' 

Uat-FNUv3t(T3team) 

Muiv*  I  8 . 0 1  5 

IF  I  ft-]  THEN  Mu;  V- 187.59 

TC  T#*.'?  TUCM  (TiC. 

▲  1  l(l^«.  iltM«  Cte.WW 

.Nfrgl-100#(Pvapl-Ps3t  >/'.  Pvapl-<  1  . Mwv/28 . 96  )  I’Paat  ) 

Nfrg2-100*<Pvap2-Psat  )/(Pvap2-( 1 .-(Muv/28.96  > )*P3at ) 

IF  Ic3v-2  then 

PRINT 

PRINT 

PRINT  USING  "10X,""Data  set  number  -  "".DO' 

OUTPUT  709; "AR  AF40  AL40  UR5" 


PRINT 

IF  Iov-2  and  Ng-1  THEN 
PRINT  USING  ''IX.'"'  Psat 


"go 


Ptran  Tirana  TaaKP) 


NG1X 


2158  PRINT  USING  "IX 


(kPa)  <kPa)  ( kPa ) 


(C) 


(C) 


Molal 


NG2S' 


Molal" 
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21B1  PRINT  USING  ” 1 X .3< 3D . DD .2X ) .2( 3D.DD .2X  )  .2( M3D.D  .2X >“ ;Pka .Pkg  ,Pkp  .Tvap2  .Taa 

t  .Mfngl  ,Mfng2 

2!B4  PRINT 

21B7!  END  IF 

2170  IF  Mfngl >.5  THEN 

2173  BEEP 

217B  IF  In-I  AND  Ng-1  THEN 
2179  BEEP 
2182  PRINT 

2195  PRINT  USING  “ 1 0X  . ‘ “Energi ze  the  vacuum  system 
2  log  BEEP 

2191  INPUT  "OK  TO  ACCEPT  THIS  RUN  ( 1-Y .0-N >?"  .Ok 
2194  IF  Ok-0  THEN 
2 1 97  BEEP 

2200  DISP  "NOTE:  THIS  DATA  SET  WILL  BE  DISCARDED!!  " 

2203  WAIT  5 

cnrn  t  Tom 

2209  END  IF 
2212  END  IF 

1  c  CMn  Tc 

^  I  O  ^  I  ^  W  A  i 

2218  IF  Ir5“l  THEN 

2221  IF  Fr.<10  OR  Fr!>100  THEN 

2224  Ifrt-0 

DCCD 

4.  ^  4.  ■  WWW! 

2230  INPUT  "INCORRECT  FM  ( 1 -ACCEPT  ,0-D£LETE{ DEFAULT  ))". I fn 
2233  IF  Ifm-0  THEN  1804 

7*?7C  CKin  TC 

41 

?77Q  CMH  TC 

b  W  W  <  4  < 

2242'  ANALYSIS  BEGINS 
2245  Ti 1 -FNTvsv58( Enf ^ 3  ) ) 

2248  Ti2-FNTvsv55{ Enf ( 5  ;  ) 

224?  Tc 1 -FNTv5v59( Enf ( 4  )  ) 

2248  Tc2-FNTv5v55^ Enf { 5  ) ) 


2249 

T i 1 .T , 1 

14  1  <41  4.  •  W  < 

1  5 

^  4.  ^ 

T  ,  7»T  , 

14  4.  14^  4.  >  4^  • 

1  c 

225  1 

To  1 -Tc 1 “273 . 

1  5 

2252 

T«7«,Tp7-.7-7*X 

1  W  4.  1^4.  4.  I  W  • 

1  5 

2254 

1«Tp1»T« 

<  W  C  4  ■  •  W  1  14 

1 

2255 

Tdcl2-To2-Ti 

2 

2256 

Tdel3-T2-T1 

?7C'7  C*  1  .C^f  ^  7  N4.C4- p  /7« 

2259  Dtde-2.5931E-2-1 .504G4E-B*Etp1  +  ! . 2 1 701 E-l 0*Etp 1" 2-5 . 1 lB4E-15*Etp1  ''3+3.2201 

E- 1 9^Etp 1  "4 

2259  Tri3-Dtde*Etp/10. 

2252  To3-Ti1+Tris 

2253  PRINT  USING  "IX,""  TIN1  T0UT1  TIN2  T0UT2  TIN3  T0UT3  DELT1  DELT2 

nci  T':;  TDTI  C  . 

^  ^  '  W  It  4  4.^ 

2254  PRINT  USING  "IX.""  (Teflon)  (Metal)  (Quartz)  """ 

2255  PRINT  USING  " 2X .  1  0(  2D.DD ,2X ) “ ! Ti 1 ,Tol ,Ti2 .To2  ,T 1  ,T2  ,Tdel 1  .Tdel2  .Tdel3  .Tri s 
2255  Er1-AB5(Ti1-T1  ) 
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7  C  '7 

C^7.ADC / T . 9-T 1  ^ 

22S6 

PRINTER  IS  1 

*^70 

IF  Erl >.5  OR  Er3> 

.5  AND  In-1  THEN 

'?'7'7  1 
^  •  • 

BEEP 

2272 

PRINT  "QCT  AND  TC 

DIFFER  BY  MORE  THAN  0.5 

2275 

OM  *  1 

'^273 

END  IF 

'^287 

TC  n|,- 1  *0  amq  Er1> 

.5  AND  In-1  THEN  1 

7S0 

2299 

IF  Ok  1-0  AND  Er3> 

.5  AND  In-1  THEN  1 

790 

9  9  Q0 

Er2-ABS ( ( T2“T 1  )-( 

Tris ) )/(T2-T1 ) 

9907 

IF  Er2>.05  AND  In 

-1  THEN 

22  9E 

BEEP 

'^''59 

pDTMT  "QCJ  AMH  T- 

PILE  DIFFER  BY  MORE  THAN 

9702 

012-1 

IF  01-2-0  AND  Er2'’ 

^  0C  AMO  T  .«»| 

1 790 

2309 

END  IF 

97  11 

PRIN'^ER  IS  701 

7  3  1  9 

tmd.jt  "COOLi^NT  TEMF .  RISE  MEi'^S .  ( 0“ 

'TEFLON 

T  ^  (Va 
r'y.  ]  A 
r  -^.1  c. 


T  C  T  +  THE^* 

C  C  T  KiT  i  IC  T  KiC.  “  t  V 


j'SING  SINGLE  TEFLON  THERMOCOUPLE  FOR  COOLANT  TEMP.  RISE" 


T  1  ^  *  T  ,  1 
T  1 

ens-  if 

TC  then 

DC  T  r,tT  I  IC  T  MC  '•  1  V 


iitTM.D  cTK-.c:c  MCTAi^  SHEATHED  thermocouple  FOR  COOLANT  TEM 


F  ■ 

^  Q 


“>  *?  T 

T  ^  ^  ^  9 

2335 

E  2  1  ^ 

-77C' 

TC  ;*r*2  Then 

^  w'  'D 

pcTKiT  ijCTMC  ■■  1  >  '■  "USING  OuARTZ 

THERMOMETER  FOR 

COOLANT  TEMP.  RISE" 

2341 

T  1  1  -T  1 

“’’^4  4 

T*^c*T2 

2347 

CMO  TC 

A  1 

2349 

TC  j  +  THEN 

9  7^G 

PRINT  USING  “ 1 X  ""USING  THERMOPILE  FOR  COOLANT 

TEMP.  RISE . 

9  7  C  0 

T  1  1  -T  :  1 

^35  1 

T  rs.T  ^7 

9  7  C 

CMO  TC 

9  7C7 

T  .  .  /  T  ^  -*.T9  ^  ^  «  C 

iQv^  'I'A-ibW'  iW 

2354 

1  f  t  -  0 

2355 

On^.  .•C^lPn^.  .r  ^ 

2359 

Rhoui-FNRhow(Tavg  ) 

2352 

Md-(E.7409^*Fn+13.027  )/1000. 

2355 

1 ^0^CC— 1 ^0CC^^C  — 

.252E-6»T1 : "2  )/l . 

0037 

9  7C  Q 

Mf-Md/Rhow 

977  1 

Uw-Mf /( PI*D: "2/4  ) 

2393 

IF  Inr,-0  AND  Uu;>.5  THEN  T2c-T2 

c-<-2.73E-4+1 .75E 

'-4  *Uw  +  9 . 35E“4*'Jw"  2- 
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Vui'-3) 

2399 

IF  Inn-1  THEN  T2o-T2o-( -B  . 44E-5+1  .71E-3*‘Ju)+4. 45E 

2407 

IF  Inn-2  THEN  T2o-T2o-( -3 . 99E-4+2 .75E-3*'v*w+1  . 45E 

2408 

IF  Inn-3  THEN  T2o-T2o-<  9.57E-5+1  .23E-3*‘Jw+1  .08E- 

2413 

Q-Md*Cpu« ( T2o-T 1 1  ) 

241S 

Qp-Q/(PI*Dc*L) 

2417! 

ITERATE  TO  FIND  OUTER  WALL  TEMPERATURES 

2418 

Tumi -Tun 

2418 

Tun2-Tun 

2420 

Tun  1 -(  Tuni +Tun2  )/2. 

2421 

Ct-Q«LOG<Do/Dr  )/( 2 . ♦PI*Kcu*L  ) 

2422 

Tun2-Tun+Ct 

ri  A')!!. 

IF  ABS<Tun2-Tun1  >>.001  THEN  GOTO  2420 

2424 ' 

ITERATE  TO  FIND  INNER  WALL  TEMPERATURES 

2425 

Tun3“Tun 

2428 

Tun4-Tun 

2427 

Tun3-(  Tun3+Tun4  )/2 . 

•7 

i.  H  ^  w 

Ct i-Q*LOG(  Dr/Di )/( 2 . *PI ^Kcu*L  ) 

24*^9 

Tun4-Tun-Ct i 

2430 

IF  ABS(Tun4-Tun3 >>.001  THEN  GOTO  2427 

243 ) 

T  unc-0 . 

2432 

T. .  ..n 

•  Mtl'l  , 

24^3 

FOR  1-0  TO  Inut-1 

2434 

Tuc( I )-Tu( I >+Ct 

2455 

Tui( I )-Tu( I >-Cti 

243B 

Tunc-Tuno  +  Tuo(  I  ) 

2437 

IWI<A  -  A  ■ 

2438 

MCyT  1 

2439 

Twr^iC^Tuinc/  Inujt 

2440 

T.  ,  *T.  T  «.  .4 

IW<<A  Al’Wk'W 

2441 

PRINT 

443 

DOtMT  nCTMC  "IV  "  "  D  ^  ,  4  <  ^ 

ll«XI4l  WWA'’*W  l\./S4VAWII 

: 

C  . . 

2444 

) 

PRINT  USING  "  iy  '‘Local  outer  wall  tenps. 

(DEG  C 

2445 

PRINT 

2447 

PRINT  USING  " 1 X “Average  outer  and  inner 

uall  ti 

T  l.tfM 

W  ,  I  w<-| 

A 

"^449 

PRINT 

2450 

Ift-0 

245 1 

Ku-FNKu(Tavg> 

2452 

Muu-FNMuu( Tavg  ) 

2453 

Muui-FNMuu(Tuni  ) 

2454 

Rei“Rhcu»Uu*Di /Muu 

2457 

Pru“FNPru( Tavg  ) 

2458 

Dtu-Q*LOG(Dc/Di  )/( 2«PI*Kcu«L  )*.5 

2459 

Lntd-(T2o-T1 1 )/L0G( ( Tuni-TI i  )/( Tuni-T2o > ) 

2480 

Per-in-PI*Di 

2481 

Surfai-PI*Di*L 

24821 

PRINT  Rei  ,Pru  ,Dtu  ,Lntd .Per in .Surfai ,Q 

rz) 

rz) 

■Z) 


eriC3 


12  3  4 

,e(DD.DD,lX)";Tuic(* 

“ "  .2(DD.DD  .IX  )"  sTuPi 
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24B7  ftnareal-P!*(Do''2-Di'2  >/4. 

24E8  ftr;area2-PI*(Dr'2-Di  "2  >/4. 

24B9  Hi-Q/(Surfai*Lmtd> 

2470  Ff  1 -<  Kcu*Perir*i*Anareal  )''.5 

2471  Ff2-(KcLi*Pcrir^*Anarea2  I"  .5 
2472'  PRINT  ftnareal  ,Anarea2  .Hi  ,Ff 1  .Ff2 

2474  Hi!-Hi 

2475  Hi2-Hi1 

247B  Alanl  •(  (  Hi2*Peri(»i )/(  Kcu*Anarra  1  ) .  5 

2477  Alam2-(  (  Hi2*PeriM  >/ <  Kcu*Anarea2  ))''.S 

2478  Ff3-FfMHi2''.5*(Twni-Tli  ) 

2479  Ff4-Ff2*Hi2''.5*(Tuimi-T2o) 

2480  F!incG“Ff3'»FNTanh(  Alanl  *L  1  ) 

2481  F'jr.c7-Ff4'*FNTanh(  Aiam2*L2  ) 

2482  Fiir!cS-Hi2*£urfai*(Twni-Tavg ) 

2483  F  jncx-FiincE  +  Func7+Func8-Q 

2484  DfuncG-(  .5*FuncE/Hi2  )+2 . *Ff 3*L1 /( 1 . +FNCoah( 2 . ♦Alanl ’L 1  )) 

2485  DfLinc7-(  .5*Func7/Hi2  )  +  2  .  *Ff  4*L2/<  1  .+FNCo3h(  2  .  *Alan2*L2  )  > 
2485  DfLir,c8-SLirfai4(Twni-Tavg) 

2497  Dfur,CA-DfuricE+Dfunc7+Df  unc9 
2488  Hil-Hi2-'.Fuinc.’,/Dfijncx  ) 

2400  7p  i^DC(  1  _Uj^2  >>  . 05  THEN  GOTO  2475 
2490 

25001  ppT|VjT  U; 

2521  Hfg-FNHf 5' Is t can  ) 

2527  TfiIr^"T5tear/3  +  Tu;no*2/3 
2530  K  f • F  NK  w ( T  f 1 1 n  ' 

2533  Phcf  "Fr'iPhc;.,''  Jf  ilr.  ) 

W|--'  . 

2535  Mu  f  *FNMljw T  f :  1  ) 

2572  01*500 

2575  Q 1  cs 5*0 1  /  ■  1  00'25  )* '  T s t caw-Trccri )  *  TO  BE  MODIFIED 

-’C'7C  +  c  \ 

25S4  Mdv*0 

7CQ'7  Dw/^  1  *  1  ^  "'O  /C  7C 

2590  H3c-Cp5C‘( Tateam-Tcon  ) 

2593  Ndvc-; ( Bp-Olcss  l-Hdv^Hac  )/Hfg 

2595  IF  ABS( ( Mdv-Mdvc  )/Hdvc  )>  .01  THEN 

200  0  N  d  V “ ( N  d  V + N  d  V  c  )  *  .  5 

2502  GOTO  2593 

2G05  END  I F 

2B09  Mdv“( Mdv+Ndvc  )•  . 5 

2E11  ‘Jg-FNUvat(  Tatear. ) 

2514  Uv-Mdv*‘Jg/Ax 
2524  Tdcf-Taat-TuP)c 

2534  Ho-Qp/Tdcf 

2535  Pr-Muf ♦Cpf /K f 

2535  Retp-(‘Jv^Dc*Rhcf  )/Muf 
2637  Nua3-(Ho»Do >/Kf 
2538  Nurel-Nuaa/Retp'' .5 
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2539  Ff-9.81*Do*Muf*Hfg/<yv''2*Kf*Tdcf ) 

2540  Hnu5s-.728*(Kf ''3*9.81*Hfg*Rhof*2/(Muf*Do*Tdcf  )  >''.25 

2542  PRINT 

2543  PRINT  USING  “IX,'"'  Uinf  Tdcf  Qfl  Ho'"'“ 

2544  PRINT  USING  “ !  X  ,2^  DD .  DO  .2X  )  .2(  M2 . 3DE  .2X  )  ”  s'.'v  .Tdc  f  ,Qp  .Ho 

2545  PRINT 

2545  PRINT  USING  "  1 X  . Nu  Retp  NuRe  Pr  F  Hnus 


2548  PRINT  USING  “ 1 X  ,2( M2 . 3DE  ,2X  )  .2( DD.DD .2X  )  ,2( MZ . 3DE  .2X >” ;Nus3 .Retp .Nurel ,Pr , 
Ff  ,Hno3  = 

2554  Rat-Mow/M'uiUii 
2554  ftr!L:s-Hi*D:/Kw 

2574  ‘Ju;*Di»Rhow  >/M'ju 

2575  ftpr5-Muw*Cpui/Kui 

2575  StcG-ftnu5/(ftrs3''  .8*ftprs".3333*Rat'-.  U) 

2577  Nupr-ftnus/  ■:  ftpra •' .  3353*Rat  '•.14) 

2070  pp  T  KJT 

2530  PRINT  USING  "  ’  X  . Ucool  Rec  Hi  Nuc  Stcoef 


2591  PRINT  USING  " I  X  .5 ( M2 . 3DE  ,2X  ) " ; Uw  .Rei  ,Hi  .ftnua  .St co  .Nupr 

2592  INPUT  "CHANGE  TCOOL  RISE'’  !“Y  2“N''  Itr 

2593  IF  Itr”l  THEN  GOTO  2312 

2594  OUTPUT  @Fi lei ; Tuo( ♦  ) 

2704  IF  Iri“1  THEN 

2707  deer 

2710  INPUT  "OK  TO  STORE  THIS  DATA  SET  (1*'Y  0"N>'^"  Oka 
2713  IF  0f.a“1  THEN 

2739  CMH  Tc 

2799  beep 

2791  TMDijT  "WILL  THERE  EE  ANOTHER  RUN  (  1  "Y  0”N>'’'‘  Igo 
2'’B4  Nrori-J 

2797  Tc  tucm  1777 

2770  ELSE 

TC  I/M,.,,,,  TUCM  I'm 


TO--  DCCD 


2935'  INPUT  "ENTER  RLOT  FILE  N,AME''  FplctS 

"'OTQl  ACCTCM  KCil-i  TO  C-1-1® 


\  cno  T.1  TH  M 


u  i  *  1  lit  Mr'' 

tv  «  I  <W 


9C/1C  I  CMTTC 


CMTTC 

w>«iwv  jtlw 


29BB •  NEXT  I 

3001  PRINT  USING  “10X,""NOTE:  “",22,""  data  runa  were  atored  in  file  "",10A";J, 
D_file$ 

3007  BEEP 

■7(^11711  TC  Tn«-'  TUCW 

3013  PRINT 

3015'  PRINT  USING  "10X,""NOTE:  "",22,""  X-Y  paira  were  atored  in  plot  data  file 
,10A";J  ,P_file$ 

3019'  END  IF 
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m/  W  I  W 

ftSSISN  ®File  TO  * 

•7ffl'7C 

W  V  I  w 

ftSSISN  eFilel  TO  * 

3079  * 

ftSSISN  @File4  TO  * 

3092 

IF  130-2  THEN  CALL 

3092 

IF  130-3  THEN  CALL 

*7  1  t  c 

END 

3119 

DEF  FNPv3t(Tc) 

■7^91 

w  1  ^  • 

COM  /Fid/  I  ft 

3124 

n  T  M  L'  f  Q  ^ 

W  ^  1  •  «  \  \  W  < 

•z  1  27 

Tp  Tf+-0 

*7  1  A  W 

DATft  -7.591234564 

3 1  33 

HATfij  4,15711732  20 

3 1 3E 

READ  K  (  ^ ) 

3 1  39 

T*{  Tc  +  273. 15 )/S47. 

3142 

Cl  1  w  0 

3145 

FOR  N*0  TO  4 

3 1  49 

C.  +  /  M  \  *(  1  -T  N 

won  woii  *  •'  V  /  •  \  »  It 

*^151 

MCVT  M 

»  *4_  /  <  »  »  • 

3154 

Br*Sofn/ (  (  1  +K (  5  )♦ 

“7  1  57 

Pp.cyp / Qp  ) 

*7  1  c  o 

■  w  <1^ 

F*22 1 20000*Pr 

" )  93 

END  IF 

1  cc 

Tp  If+-1  THEN 

“7  J  CQ 

Tf*Tc»1 .9+32+459.5 

'7179 

P»'0'‘<  33.0555-4330 

7  1  75 

p»p*t0)37c/ 14^595 

7  1  70 

END  IF 

7  1  w  1 

T p  ] f  *  »2  Then 

7  1  Cf  /I 

A*°. 394595-3055 . 1  / 

7  1  07 

p«  !  *3 . 32 *  1  0  '■  A 

7  1  af?i 

•o  1  w  V 

END  IF 

3 1  93 

pCTlipiVj  p 

7  1  QC 

•  kJ  w 

FNEND 

3 1  99 

nCC  T  \ 

woi  .  ..11  n  y  \  / 

320*^ 

CON  /Fid/  Ift 

3205 

IF  Ift-0  THEN 

3209 

Hfg-2477200-2450»( 

77 1  1 

w  ^  ■  ■ 

CMR  TC 

0  1  «  W  *  t 

77  1  ^ 

IF  Ift-1  THEN 

7  7  \  7 

T  f *T  *  1 . 9+32 

"^220 

Hfg-7.0557957E+1-T 

7777 

U«^«U«^»777C 

n  <  ^  1  >  1  ^  4.  w  k  w 

3225 

END  IF 

3229 

IF  Ift-2  THEN 

3232 

TK-T+273. 15 

3235 

Hfg-l .352E4E+6-TK* 

3239 

END  IF 

3241 

RETURN  Hfg 

3244 

FNEND 

3247 

DEF  FNMiJw(T) 
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COM  /Fid/  Ift 
IF  Ift-0  THEN 
Fi-247.8/(T+133.15) 

Mu-2.4E-5*)0*ft 
END  IF 

IF  Ift-1  THEN 

Mu-8.9629819E-4-T*( 1 . 1 094B09E-5-T*B . SBB929E-9  > 

END  IF 

IF  Ift-2  THEN 
TK-!/(T+273.15) 

Mu-EXP' -1 1 .0’79+Tk*( 1 . 744E+3-Tk*( 2 . 9033SE+B-Tk* 1 . 12B61E+8  >  >  > 
END  IF 
RETURN  Hit 
FNEND 

nCC  CMUw,* / T+ \ 

COM  /Fid/  If*’ 

IF  Ift-0  THEN 

D.CKiDw.  i  /  T|  \ 

T-Tt  +273 . ] 5 
X  - 1 5  0  0  /  T 

pi -1 /( 1 +T. 1 ^C_4 ) 

F2-(  1  “EXP(  -X  )  )''2 . 5*EXP<  X  }/X'^ .  E 
9- .001 5*F1 - . 000942 *F2- .0004SS2*X 
u  .2<‘P/(  46  1  ,  52*T  ) 

Urn  (  1  J.  /  1  4.  2  4  g  *  K  )  '* .  5  )  /  K 


IF  Ift*l  THEN 

Tf*T+*1 ,0+32 

v‘*13.9553E7“Tf*(  .  161272E2”Tf*^.  1729''  90E“4  ) 

I  l«i  1  /  1  C  <7)0 


IF  Ift-2  THEN 
tl  mT*  i  c 


i .  9E*Tf'  /P 


RETURN  U 
FNEND 

n£f  PMQ"-/  T  ) 

rnM  /c  ^  M  /  r  f  * 

ww'l>  /  I  AW/  Al  V 

TT  +  TUC^^ 

Al  AIV 

Cpw-4 .2 1 1 Z0959-T#( 2 .2682BE-3-T*( 4 .423B1E-5+2.71429E-7*T  > ) 

END  IF 

IF  Ift-1  THEN 

Cpw-9.2507273E-1 +T*( 8.3400433E-4+! . 7207792E-B*T ) 

END  IF 

IF  Ift-2  THEN 
Tk-T+273. 15 

Cpu;-4. 18B8*(  1  . 5984E-2  +  Tk (  3 . 35083E-3-Tk*(  7 . 224E -B-Tk*7 . 6 1 749E-9  )  )  > 
END  IF 
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lAi  c>j 


f*9.2494SE^4-  +  TK^(E.2795E“7  +  Tk^(9.2444E”10-i'Tk'’»3.0E7E“12)) 


RETURN  Cpw^!000 

“^'403 

FNEND 

340W 

DEF  FNRhcw(T) 

3409 

POM  /Fid/  Ift 

34  ?  2 

IF  lft"0  THEN 

■741c 

RO-999.52946+T*' 

7410 

END  IF 

■547  1 

IF  Ift*l  THEN 

*7  474 

Pq.i  E207479E+3- 

3427 

CKin  Tp 

"^430 

IF  I f THEN 

343*^’ 

Tl.-«Tx7'73_  lE-339 

343E 

lifvO  74O4QC-4XTI 

7  4  7  Q 

P  D  ■  1  ■  U  f 

*^44 

END  IF 

prTtJDM 

"^44^ 

FNEND 

34E  I 

nCT  T  s 

•  •«<  I  WW'  ‘  • 

"^4^4 

T  \*pMi 

3457 

RETURN  Rru; 

34£0 

FNEN*^' 

74C  ' 

nrp  j  \ 

74c  c 

PPM  Tf  + 

74c  5 

IF  I  f  t  *0  THEN 

7477 

X  •  (  T-^273  .  >  5  )  /  27 

7475 

L  ^  922477y  •  ■  2 

“?  4  7C 

c  Kj  n  T  r 

4  ' 

TC  Tf*»1  TuZt\ 

7  4  c  4 

^  . 209523SE-2 

7  4  c  7 

END  IF 

"^490 

Ift«Z  Then 

34  93 

T 1  « T '*■2 73  .  ■  5 

34SE 

L  .  .W.A  1  OC  OC-/ •  /  C 

<^Wk>  **f-WVJWW  V  ^W 

3499 

CKir.  TC 

^  I  <  W  A  > 

3502 

RETURN  Ku; 

3505 

CMCMR 

•  •  1  «  by 

750s 

DEF  FNTflph^  X ) 

35  •  1 

D  «  C  y  0  (  V  1 

7C  ’  4 

n  »  E  X  P  C  ”  X  ; 

7-17 

T.r>k»^c-.n 

1  w  >  •  •  >  <  '  wr  '  '  ur 

7C  »7 

PE"^URN  TapH 

3523 

CMC  Ki  n 

^S25 

n  C  p  p  K1 T  M  J  \ 

7C7  a 

CON  ^ C  c /  C  (  7 ) 

T«r  (  o\  \ 

I  ^  / 

FOR  I ■ !  TO  7 

TmT^Cf  T  ^♦U''T 

•  ■  ’  W  '  A  /  w  A 

NEXT  I 

T-T-"4.733eBE-3+T*(7.S92S34E-3-T*8.e77927E-5> 
RETURN  T 
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(J  IjJ  tA)  1>J  1>J  tjJ  1>J  t)J  1>J  t>l  l.xj  1>J  l>)  (>1  t>)  UJ  t>J  l,.J  1>J  t.-l  l-l  t,-l  t>)  t>)  l.-l  W  t>4  1>J  1>1  t>i  1>J  1>I  t>a  OJ  1>J  t>l  1.>J  IaI  t>J  (>I  (.>J  W  lAl  IaI  CaI  lAl  CaI 

ID  to  ID  lO  tD  ID  tD  LD  CD  CD  CD  CD  CO  CD  CD  CO  CD  CO  CO  CO  CO  i  IJ  CO  CO  CD  CO  CD  CD  CO  CD  CO  CD  CD  CO  CO  CO  CO  CO  CO  CO  CD  -J  -J  --J  -J  -J  --J  -O  -J  -O 

CaI  r-j  —  ca  C3  ca  (9  ca  C9  e3  c3  CD  lo  co  cd  cd  co  CO  CD  CD  CD  lo  CD  CO  CO  CO  CO  CO  CD  CO  CO  -j  --J  -o  cn  tn  C.-J  ro  --  c3  CD  CD  -o  -j  -g  -g  -g  cn  cn 

-g  -g  -g  co  co  -g  cn  eg  ro  ts  cd  cd  -g  cn  tn  i.-j  to  -  ra  id  id  -g  cn  tn  c.-i  r o  —  ea  cd  cd  -g  -g  -g  -g  -g  -g  -g  -g  -g  -g  -g  cn  m  ■!:»  eg  cd  co 


RETURN  T 
FNEND 

DEF  FNTvsvESCJ) 

roim  /r,.c;p/  Tt:o/ c  c 
TmTCOf  Ok  \ 

I  •  ^  W  V  w  / 

FOR  I-l  TO  5 

T.TxTCQ/  T  \^H''  T 

NEXT  I 
RETURN  T 
FNEND 

DEF  FNCc5h(X) 
P"EXP( X  ) 

Q-EXP(-X  ) 

Cc5 h* , 5^ ( P+Q  ) 
RETURN  Cosh 
FNEND 


CON  /Fid/  Ift 
DIN  Errf  (  1  3  ) 


INPUT  "ENTER  FILE  N/^NE"  Fiie$ 

ACCTCM  «C  -  1  •  TO  C  <  ^ 

BEEP 

INPUT  "ENTER  PRESSURE  CONDITION  ( 0*U ,  Ipe 

INPUT  "ENTER  NUNBER  OF  RUNS”  Nrurt 

OOTMTCO  TC  nO] 

>  AW  iwl 

PRINT 

DO  T  MT 


1  ••  {  C  ^  ^ 


TC  TUCM 


PRINT  USING  ”  I  0X  "Pressure  Ccrditior:  {^tnospher  :c 


CMTCD 

k.l«lk,IC  \;i  AAW«AI^,AIMI 


Cc.  Ch 


in 


Tout 


ENTER  §Fils; Inwt  .Ff 
PRINT 

PRINT  USING  "10X."''Dat3 
FOR  I-i  TO  Nrun 

ENTER  @File;Bvcl  .Sarp  .Ptran  ,Etp  ,Enf(*  )  ,Fr!,Tl  ,72  ,Pv 
PRINT  Evcl  .Bar.p  , Ptran  .Etp  ,Ef<\f  (  *  )  ,Fr.  ,T1  ,T2  ,Pvap1 
NEXT  I 


ftSSIGN  @File  TO  • 

SUBEND 
sue  Purg 
BEEP 

INPUT  "ENTER  FILE  NAME  TO  BE  DELETED" .Fi le$ 
PURSE  Files 
GOTO  3909 
SUBEND 
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